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ABSTRACT 


The effectiveness of spectral analysis in studying the 
responses of single motor units was investigated for motor units in the 
iliofibularis muscle of the toad Xenopus laevts. 

Isolated motor axons were electrically stimulated with 
narrow, supramaximal pulses. Motor units which displayed a twitch-like 
contraction in response to a single stimulus pulse were classified as 
fast motor units, and those which did not display such a contraction 
were classified as slow motor units. On the basis of their isometric 
tension responses to stimulation with trains of evenly-spaced pulses, 
the fast motor units were further subdivided into three types, and the 
slow units were subdivided into two types. 

The histochemistry of muscle fibres from several isolated 
motor units was studied. The results of such studies support the 
hypothesis that each unit contains only a single type of muscle fibre. 
Each motor unit type was matched with one of five muscle fibre types. 

Spectral analysis me thods were used to study the responses 
of each type of motor unit to trains of randomly-occurring pulses. 
Transfer function and coherence estimates which related the low freq- 
eet scone oan s of the random stimulus train to the resulting isometric 
muscle tension were calculated. These estimates indicated that for 
frequencies in the range 0.5-25 Hz, the response of a slow motor unit 
to random stimulation is poorly represented by the response of a linear, 
noise-free system, and the response of a fast motor unit is similar to 
that of a second-order low-pass filter with an added time delay. 


Curves representing the responses of such linear systems were fitted to 
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transfer function estimates for fast motor units, and confidence inter- 
vals for the parameters defining these fitted curves were estimated. 
For a given motor unit, these parameters depended systematically upon 
such factors as muscle length and mean stimulus rate. Computer 
simulations of fast motor units, based on the experimentally-observed 
linear models, were used to demonstrate the effects of various modif- 
ications to the spectral analysis procedure. 
The linear models for the responses of fast motor units to 

random stimulation were only approximate and were subject to many 
restrictions. However, such models were found to be useful as a 
standard basis for describing motor unit responses. They also provided 
information about motor unit function which supplemented the results of 


studies involving stimulation with trains of regularly-occurring pulses. 
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CHAPTER uf 
INTRODUCTION 

Studies of input-output relationships for the elements of a 
control system often involve the development of mathematical models 
from experimental observations. Linear relationships can be modeled 
in this way by fitting transfer functions to frequency response plots 
based on experimental data [1,2,3], and, in many situations, nonlinear 
pelaidionahies may be approximated by such linear transfer functions. 
Attempts have been made to use this approach in studying various 
elements of the system which controls movement in animals (the motor 
control system), but the use of linear transfer functions in such 
situations is complicated by the fact that many of the signals which 
are encountered appear as trains of identical, narrow spikes [4], 
eae the behavior of a linear system is generally described in 
terms of its transfer function for continuous sinusoidal signals. 

A skeletal muscle is composed of many muscle fibres arranged 
into individually-controllable groups called motor units. These fibres 
contract when the nerve serving the muscle is electrically stimulated 
with a train of narrow pulses. A mathematical model describing the 
relationship between the electrical stimulus and the resulting muscle 
tension for single motor units would be a useful component of any 
analysis of the entire motor control system; but, if such a model is 
to be determined, a means of analytically relating the pulse train to 
the continuously changing muscle tension must be devised. Partridge 
[5], who studied the frequency response of the whole triceps surae 
muscle in the cat, approached this problem by using pulse trains for 


which the interpulse intervals were sinusoidally modulated. Frequency 
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response plots expressing the relationship between the modulating 
signal and the muscle tension were produced, and transfer functions 
were fitted to these plots. A similar approach was used by Rosenthal 
et al. [6] for the cat gastrocnemius muscle. 

Pulse trains with sinusoidally modulated interpulse intervals 
have certain disadvantages when used as inputs to measure the frequency 
response of nerve-muscle preparations. The production of a frequency 
SE plot requires measurement of the responses to a large number 
of different sinusoidally modulated signals, and this means that the 
preparation must be subjected to long periods of stimulation which 
could cause muscle fatigue. In addition, because such pulse trains 
contain power only at discrete frequencies [7], any nonlinearities 
arising from the effects of intermodulation could not be detected. 

In order to overcome these disadvantages, spectral analysis 
methods have been used to relate spike trains to continuous signals for 
both sensory [8-10] and motor [11,12]. elements of the neuromuscular 
system. For the nerve-muscle preparation, a train of randomly occur- 
ing pulses is used as the input signal, and the resulting transfer 
function estimate relates the frequency components of this pulse train 
to the corresponding frequency components in the continuous signal 
representing muscle tension. The input pulse train can be tailored so 
that it contains power at all frequencies within the band of interest, 
thereby insuring that any potential intermodulation effects will occur. 
If the system is not perfectly linear and noise-free, the transfer 
function estimate will represent the response of the linear, noise- 


free system which describes with minimum mean-squared error [13] the 
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actual system response. The degree to which this linear system 
contributes to the actual system response is apparent from the coherence 
estimate which is readily calculated from the same data used in deter- 
mining the transfer function estimate [14]. 

Several approaches are available for the spectral analysis of 
random pulse trains: the pulse train may be viewed as a realization of 
a stochastic point process and analyzed accordingly [15-18], it may be 
viewed as a train of unit impulses [19], or it may be filtered with a 
low-pass filter to produce a continuous, band-limited signal [8,9]. 

The filtering approach is attractive for the study of nerve-muscle 
preparations because special hardware to measure the occurrence times 
of the impulses is not required, and because the continuous signal 
appearing at the filter output may be processed with the fast Fourier 
transform [20-22] in the same manner as the continuous signal repre- 
senting the system output [14,23,24]. 

In considering the application of spectral analysis to the 
study of input-output relationships for single motor units, several 
questions may be posed. 

LS What spectral analysis procedures should be used for this 

application? 

2. Can the results of such studies be readily interpreted? 

3. How are the transfer function estimates related to the 

responses resulting from eecietien with trains of 


evenly-spaced pulses? 
4, How does the structure of the muscle affect the 


observed input-output relationships? 
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5. How does a study of. the input-output relationships for 
single motor units contribute to our knowledge of motor 
unit organization? 

Beginning with reviews of some of the relevant areas of muscle 
physiology and spectral analysis in Chapters 2 and 3 respectively, this 
thesis is an attempt to answer the above questions for motor units in a 
muscle of the hind limb of the amphibian Xenopus laevis. When this 
project was begun, it was not known whether spectral analysis methods 
could be successfully used in such an application, but, since that time 
the results of a few studies demonstrating the feasibility of such an 
approach have appeared [11,12,25]. Therefore, although the problems 
involved in implementing the technique are considered here, much of the 
axendedd is placed on the interpretation of the resulting transfer 
function estimates with the aim of determining whether these results 
provide enough information to justify the considerable effort required 


to obtain them. 
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CHAPTER 2 
MUSCLE FIBRES AND MOTOR UNITS 
Bigal Vertebrate Skeletal Muscle 
2.1.1 Structure 

According to current phyadieldar textbooks [27-29], vertebrate 
skeletal muscle is composed of an aggregation of cylindrical muscle 
fibres held together with connective tissue. These fibres are 
arranged longitudinally, and they often run the full length of the 
muscle. Unfortunately, most of these sources imply that vertebrate 
skeletal muscle is composed of identical muscle fibres when, in fact, 
several types of muscle fibre have been described, both in mammals 
[30-33,49] and in amphibia [34-37]. The contractile properties of 
individual muscle fibres can vary greatly from fibre to fibre in the 
same muscle, and the distribution of fibre types varies from muscle to 
muscle in the same individual. However, the structure of a particular 
muscle is remarkably similar in different individuals of the same 
species. 

Muscle fibres are composed of bundles of cylindrical fibrils, 
and each fibril is composed of many longitudinally arranged filaments 
which are organized into distinct bands [38]. A system of nomenclature 
to describe this band-like structure has arisen (Figure 2.1). The A 
(anisotropic) band is composed of thick parallel filaments 100-120 A 
in diameter and 1.6jymlong. The I (isotropic) band filaments, which 
are generally attached at their mid-point to a structure known as the 
Z line, are thinner (60-70 D in diameter). The two sets of filaments 
overlap within much of the A band, and they are generally arranged so 


that each A filament is surrounded by six I filaments; but other 
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alignments do occasionally occur. The region within the A band where 
the two sets of filaments do not overlap is known as the H (Hensen's) 
band, and the portion of the fibril delineated by two Z lines is called 
a sarcomere. This arrangement is repeated along the length of the 
ribet 
An important component of every muscle fibre is called the 

sarcoplasmic reticulum [39]. This structure is a continuous, membrane- 
limited system of tubules which form a close network around each fibril. 
The finger-like tubules which surround the myofibrils are longitudinally 
oriented, and they periodically converge to form dilated circumferen- 
tial structures known as terminal cisternae. Another system of tubules 
within the muscle fibre is in close contact with, but separate from, 
the terminal cisternae. These tubules, which appear.to be or to the 
extracellular space at the cell membrane [40], run transversely through 
most of the muscle fibre and are therefore called transverse or 
T-tubules. The arrangement of the sarcoplasmic system and the T-tubules 
is such that certain characteristic features known as triads occur in 
every pe A triad consists of a T-tubule and two terminal: cis- 
ternae arranged so that a longitudinal section through the muscle would 
reveal the T-tubule between the two terminal cisternae aligned along an 
axis aired to the myofibrils. In most mammals, the triad is 
centered near the A-I junction, but in frogs, it coincides with the 2 
line [41]. 

| As well as the more or less regularly arranged myofibrils and 
membrane systems, the muscle fibre contains many irregularly placed 


bodies. Scattered throughout each muscle fibre are structures known 
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as mitochondria which are involved in the synthesis of adenosine 
triphosphate (ATP). Glycogen and droplets containing fat act as 
sources of energy for this process. Every muscle fibre also contains 


many peripherally-placed nuclei [38]. 


2.1.2 Muscle Contraction 

A motor nerve cell (motoneuron) transmits information from 
the central nervous system to the skeletal muscles by means of action 
potentials (temporary depolarizations of the nerve cell membrane) [49]. 
These action potentials travel along a thread-like extension of the cell 
body known as the nerve axon. A structure called a synapse is located 
at the junction of a motor nerve axon and a muscle fibre [45]. When an 
action potential arrives at the synapse, a small quantity of the 
chemical acetylcholine is released. The acetylcholine causes a local 
depolarization of the muscle fibre's surface membrane, and, in most 
muscle fibres, this depolarization develops into an action potential 
which is transmitted along the fibre. A release of calcium ions by the 
sarcoplasmic reticulum is associated with the depolarization of the 
surface membrane [46], and these calcium ions catalyse a chemical 
reaction involving ATP which causes the interdigitated A and I fila- 
ments to slide together producing tension [41]. After a short time, 
the calcium ions are sequestered by the sarcoplasmic reticulum, and the 
muscle tension drops back to its initial level [47]. Thus a single 
depolarizing pulse applied to the motoneuron causes a single muscle 
twitch. If the motoneuron is stimulated by a series of pulses, the 


individual twitches tend to sum; and if the frequency of stimulation 
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is high enough, a smooth sustained muscular contraction results. This 
type of contraction is called a tetanic contraction. 

Although the above description of the mechanism for muscle 
contraction holds for a large proportion of muscle fibres, the response 
of certain amphibian muscle fibres does not follow this pattern. The 
local depolarization induced by the release of acetylcholine at the 
synapse does not result in the generation of a propagated action 
potential [48]. Therefore the muscle fibre contracts only in the 
region of the synapse rather than over its whole length, and it does 
not produce a twitch in response to a single stimulating pulse applied 
to the motoneuron, If the motoneuron is stimulated by a series of 
pulses, however, a smooth sustained contraction is observed [50]. 

The tension developed by muscles is used mechanically in 
several ways. The muscle may shorten and produce movement, it may 
oppose a constant force (such as gravity) without changing length, or 
it may oppose some external force or antagonistic muscle while it is 
being stretched. Although the muscle can develop tension over a wide 
range of lengths and external forces, normally two standard modes of 
contraction--isometric and isotonic--are used in the study of the gross 
contractile properties of muscle [52-54]. Isometric contraction refers 
to contraction at a fixed length, and isotonic contraction refers to 
contraction with a constant load. 

The human body contains approximately 2.5 x 10° muscle fibres 
but only about 4 x 10° motoneurons [55], and similar ratios of muscle 
fibres to motoneurons are apparent in all vertebrates. Therefore, one 
motoneuron must innervate many different muscle fibres. Since the 


control of a muscle is achieved via the motoneurons, the smallest 
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discrete unit which can be controlled by the central nervous system 
consists of all of the muscle fibres innervated by the branches of a 
single motoneuron. This group of muscle fibres, together with its 
motoneuron is known as a motor unit. Some small motor units contain 
only a few muscle fibres [56], but larger motor units containing 


several hundred muscle fibres are not uncommon [57]. 


Lae Amphibian Skeletal Muscle 
2.2.1 Response to Stimulation 

Differences between the muscle fibres in amphibia were first 
noted in 1944 by Tasaki and Mizutani [58]. They observed differences 
in the contractile properties of the muscle when various isolated 
motor axons were stimulated. Kuffler et al. [50,59] found that stim- 
ulation of the larger diameter motoneurons resulted in the familiar 
propagated action potential and resulting twitch, but selective stim- 
ulation of the smaller diameter motoneurons resulted in a slow local 
contraction and a non-propagated depolarization in the area of the 
local contraction. The two different responses were observed in the 
gastrocnemius, peroneus, and tibial interossei muscles, but not in the 
sartorius or adductor longus. The so-called "small nerve system" was 
especially prevalent in the iliofibularis and semitendinosus, and a 
bunching of the small nerve activated muscle fibres was noted in the 
iliofibularis. This bundle, which had earlier been called the "tonus 
bundle" by Sommerkamp [60] was composed of reddish-colored fibres. 
When the tonus bundle was subjected to the same stimulation procedure 
previously applied to the whole muscle, both propagated large nerve 


contractures and local small nerve contractures were observed. 
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Kuffler [61] continued to study the properties of muscle 
fibres innervated by small-diameter motoneurons. He found that the 
muscle fibres possessed many Seal contacts with motoneurons, that 
these contacts were distributed along the whole fibre length, and that 
the contacts on an individual fibre did not all originate from the 
same axon. He also noted that the teh ete potential was smaller than 
that of ene fast, twitch producing fibres and that propagated action 
potentials never occurred in the slow fibres. These results were 
later confirmed by Burke and Ginsborg [48]. Kuffler also concluded 
that. no individual muscle fibre can exhibit both "fast" and "slow" 
types of response. Kuffler and Williams [62] noted that although the 
slow fibres were incapable of developing large tensions, they could, 
when continuously stimulated, maintain large tensions originally pro- 
duced by twitches. 

Nasledov [63] studied the electrical excitability of fibres 
in the sartorius and in the tonic bundle of the iliofibularis. He 
found that the sartorius was composed exclusively of fast fibres and 
that the tonic bundle of the iliofibularis contained two types of fast 
fibres in addition to the slow fibres. One of these fast fibre types 
had properties similar to those noted in the sartorius, but the other 
type of fast fibre did not always develop an action potential in 
response to a single stimulation of the nerve. This type of fibre 
probably corresponds to those noted by Ralston and Libet [64] and by 
Kuffler and Williams [62] in which the development of action poten- 
tials is facilitated by stretch or multiple stimuli. Orkand [65] 


proposed that the reason for the low excitability of such fibres is 
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the small quantity of transmitter released at the motor endplate by an 


individual nerve impulse, 


2eeue JOLLUCtUre 

Although much of the early evidence for the existence of more 
than one type of muscle fibre was obtained from studies of a muscle's 
response to electrical stimulation, observations by microscopists also 
indicated the existence of at least two structurally different types of 
muscle fibre. Kruger's [66] classification of amphibian muscle fibre 
structures into two types, frtbrillenstruktur and felderstruktur, was 
confirmed by Hess in 1960 [67]. Ftbrillenstruktur fibres contain lipid 
droplets, have many sarcoplasmic granules and possess individual fibres 
surrounded by sarcoplasm. The Z lines of such fibres appear straight 
across the width of the fibrils. Felderstruktur fibres, on the other 
hand, contain no lipid droplets and fewer sarcoplasmic granules. Their 
fibres appear to clump together, and the Z lines appear wavy across the 
width of the fibre. 

Two distinct types of motor end-plate (junction of nerve and 
muscle) were observed by Gray [68] in 1956 and this observation was 
confirmed and expanded upon by Hess [67]. Hess stated that the nerve 
endings were either of en plaque or en grappe type. En plaque endings 
are branching, ribbon-like structures originating from large-diameter 
axons, whereas en grappe endings appear as clusters of tiny balls 
originating from branches of the smaller axons. Variations between the 
en plaque endings of fibres within the tonus bundle and those not con- 
tained within the tonus bundle were noted. kn grappe endings were 


found only on fibres within the tonus bundle,and many such endings 
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were observed closely spaced on each muscle fibre. Generally only one 
en plaque ending could be seen per muscle fibre, and when more than one 
such ending was observed on a single fibre, the multiple endings were 
widely spaced. Both types of motor ae could not be seen on the 
same fibre. Direct observation showed that muscle fibres with en 
plaque endings are of ftbrillenstruktur type only, and that en grappe 
endings occur on fibres of felderstruktur type. It thus appeared that 
the fast muscle fibres are of fibrillenstruktur type with en plaque 
endings and the slow muscle fibres are of felderstruktur type with 

en grappe endings. 

Further structural differences between slow and fast muscle 
fibres were documented by Peachey and Huxley [69]. They observed that 
the fast fibres are composed of myofibrils well-delineated by sarco- 
plasmic elements. The M band (a swelling of the thick filaments 
occurring in the middle of the A band) is well-defined, and triads are 
visible. In contrast, slow fibres are composed of larger ribbon-like 
myofibrils, often fused together; they have no visible M bands and no 
triads; and the Z bands of such fibres are thicker. In addition, the 
slow muscle fibres possess more mitochondria than do the fast fibres. 
Page [70] agreed with most of these observations, but, unlike Peachey 
and Huxley, she was able to detect triads occurring at every fifth or 
sixth Z line in the slow fibres. 

Smith and Ovalle [34] have carried the structural classifica- 
tion of muscle fibre types a step further by staining serial sections 
of the same muscle for two different enzymes, succinic dehydrogenase 
and myosin adenosine triphosphatase, and correlating the histochemical 


profiles with the fibres’ ultrastructural characteristics. With this 
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technique, they were able to distinguish five types of muscle fibre 

that they have designated numerically. A summary of their results is 
presented in Table 2.1. The sartorius muscle, which contains no slow 
muscle fibres is composed solely of types 1 and 2 fibres as is the outer 
portion of the iliofibularis. The centre portion of the iliofibularis 
is composed of types 3, 4 and 5 muscle fibres. Smith and Ovalle feel 
that the muscle fibre types designated 1, 2 and 3 belong to the 
population of fast muscle fibres and that types 4 and 5 are slow muscle 


fibres. 


2.2.3 .Relationship Between Structure and Function 

The functional properties of the various types of muscle 
fibres in the iliofibularis were studied by Lannergren and Smith [35]. 
They found that the fibres outside the tonus bundle possess a moderate 
content of lipids and oxidative enzymes, and that these fibres fatigue 
rapidly if stimulation is maintained. They also found that some of the 
fast fibres within the tonus bundle possess a high content of lipids and 
oxidative enzymes and do not fatigue rapidly. The oxidative enzyme and 
lipid content of the slow fibres was found to be very low. Similar 
histochemical results were obtained by Asmussen and Kiessling [71]. 
They identified two types of fast muscle fibre and observed that the 
tonus bundle of the iliofibularis is composed of two fibre types--— 
slow ftaree with an extremely low content of oxidative mitochondrial 
enzymes and thin, fast fibres with many mitochondria. 

Recently, Zagorska [37] has classified fibres in the eye 
muscles of various species of frogs and toads into four groups. This 


classification, based on structural criteria, includes three fast 
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muscle fibre types and a single slow fibre type. Species-related 
differences were noted in the proportions of the various fibre types in 
each eye muscle, and these differences agree with the proposition that 
the muscles are constructed to suit the particular way of life of each 
species. 

A study by Smith and.Lannergren [72] of the properties of 
individual Onay units in the iliofibularis identified four types of 
motor units--two fast and two slow. The fast motor units composed of 
fibres in the outer part of the muscle fatigue rapidly in comparison to 
the fast motor units near the centre of the muscle. The two slow motor 
unit types can be differentiated according to the rate of rise of 
tension in response to a train of stimulus pulses. Both types of slow 
motor unit are located near the centre of the muscle. The fact that the 
location of single motor units within the muscle corresponds to the 
location of a given type of muscle fibre indicates that motor units may 
possibly be homogeneous in composition. Further support for this pro- 
position is provided by the fact that the functional properties of 


single motor units and individual muscle fibres are remarkably similar, 


2.3 Summary 

Skeletal muscles in vertebrates are composed of muscle fibres 
which are not all similar, either in terms of their structure or their 
function. Muscle fibres in frogs and toads can be broadly classified 
into two groups--fast and slow--and evidence indicating the existence of 
a nets of five sub-classifications has been presented. The muscle 
fibres are organized into motor units having a variety of functional 


properties. Because each motor unit appears to be homogeneous in 


16 


is) aebi nor bid. Ja! (5%) usrqyoapkd: bas ime vd 


lo 29999 quot ertidwsbt pisabudlto2ts eds nf esting x0 cme 


ty Lgeoemoo ‘etlan rolJom Jee 2) wits ows bas yan? owl--eshnw 
5 foo ads 3 tm IJiBg Tague erty nit soya 

sQ76e waln owt SO? .slseym) oa 20 Srintso Sild  tKS cs Soson a weld “ 

te, seis Ym 3961-00097 -gakbacone Lovarepowsttin wd ies cone jing | 

dole tn peqys dj0e “Jeeetug eutomise io tet? & OF sénoqam ak notens3 . 

ofa S0df fo53 si? aise silt jo 9778989 tad y4snm, busnzol ak ealiaacea 

diff od dbnoqesi 109 sive 414, Aluiw ering wc F tgate i! oe 

vem etine 35j3a8 JS0di PoI56 slbat agi ri $foaum To Savi revtg i 46 nataapel - 

aty ehrd xOd PIOqque santret nubateoqaas: at syosnqgoaed af —tdtateg! : 


64 copdvianiios ni elbiqst aug tsG3 


26 wabiseqoi7 Lsantsonet sry ons 1980 sty yd bebivesg ak 


aeLlats: vi fidigees sy) ox: eegdls sToemn Lavbivihnt bk esta s920m 


. 
eaidt? olosum to bscoquo2 are esdardsirey i" o sa 
var) 1 sauIsuxte 4Is\t1 Io. sa1s2 at aia wat tabs 1 ioe a Ts | pe 
beljtsesis uihsosd au; fies, Bhaot i esvvt?, aloe: : 
( oe 
ty somezatas: oly. annie ea igen 2 ; 
star at boxnaeoxy oaod i arse te tt 


composition, the possibility of a one-to-one correspondence between 
motor unit types and muscle fibre types is evident. The demonstration 
of such a relationship would be a worthwhile expansion of the current 


knowledge concerning the functioning of amphibian muscle. 
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CHAPTER 3 
SPECTRAL ESTIMATION 


3.1 The Power Spectrum 


The autocovariance function RT) of a stationary random 


process x(t) is defined as 


T/2 
RE = ae tf x(t) x(errar (Se15 
-T/2 


A random process is stationary if its statistical properties are 
unaffected by any translation of the time origin. The Fourier trans- 


form of R(T) is known as the power spectrum of x(t): 


es = ecole eae dt (3.2) 


This quantity is a measure of the power per unit bandwidth (in Hz) 
centred at frequency f. Although the Fourier transform of a signal is 
generally a complex quantity, the power spectrum is real because Ry (7) 
is symmetrical about the T=0 axis. 

In order to evaluate the power spectrum of a signal with 
Equation (3.2), an infinitely long sample of the signal is required. 
In practical situations, only samples of finite length are available, 


and the power spectrum must be estimated from such samples. 


3.1.1 Power Spectrum Estimation 
In 1958, Blackman and Tukey [73] described an indirect pro- 
cedure for calculating spectral estimates based on the Fourier 


transformation of autocovariance functions evaluated over the sample 
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length, This procedure was cans until the development of the fast 
Fourier transform algorithm in the middle 1960's [20-22,92]. With this 
algorithm, spectral estimates can be calculated directly by Fourier 
transforming the finite-length sample of the process. Because it is 
much faster than the indirect procedure, the direct procedure is now 
used almost exclusively. 

A finite sample of the random process x(t) will be denoted as 


Xp (t) where 


- hxCEy Ostet 
Xp (t) = \ 0 <0, t>T | (Si a3) 


A direct estimate of the power spectrum of x(t) based on Xp (t) is 
a it 2 : 
SCD ae Gay (3.4) 


where x, (£) is the Fourier transform of X(t) F239 hy 

If S x (FD is approximately constant over a frequency band of 
width 1/T, the expected value of S x(t) is Sex (FD [23]. However, as T 
becomes large, the variance of the spectral estimator approaches Baked 
rather than zero [23]. The probability distribution of the spectral 
estimator S (ED is proportional to that of a chi-squared distribution 
with two degrees of freedom [23]. 

Because the variance of Bt) does not approach zero as the 
_ record length increases, S tf) is an inconsistent estimator of the true 
power spectrum, and often the variance of such an estimate is so high 
that the estimate is useless. A further refinement of the analysis 
procedure is required in order to eliminate this inconsistency and 
reduce the variance to an acceptable value. If n independent 
evaluations of Sf) are calculated and then averaged, the result is 


called a smoothed spectral estimate. These n independent evaluations 
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may be obtained from X(t) by dividing this record into n non- 
overlapping sub-records each of length L and then calculating a 
separate spectral estimate for each sub-record. For such a situation, 


n=T/L and the smoothed spectral estimate Booker) [74] is 


bee 


S (£) (3.5) 


xX 


Sil 


He 
I 


Sey (F) = 


p= PENS 


A particular sub-record x, (t) may be viewed as the product of 


x(t) and Wy (t) where 


fl eystst> 
Wy (t) ‘5 ou (6) 


and to-t, =L. Wp (t) is known as a rectangular data window, and its 
Fourier transform, Wp (£) is called a frequency window. Multiplication 
of x(t) by a data window is equivalent to the convolution of X(f£) with 


the corresponding frequency window [75]. That is, 


co 


K(f) =_/ W, (£-g)-X(g)dg (3752) 


This convolution results in a leakage between frequency bands which is 
generally undesirable, so the shape of the data window is often modified 
to reduce the leakage. If the abrupt jumps at the ends of the rectan- 
gular data window are smoothed out, the side lobes of the resulting 
frequency window will be lower than those of the frequency window 
corresponding to the rectangular data window, and the leakage between 
widely separated frequency bands will therefore be less [23]. Some 
commonly used data windows and their corresponding frequency windows are 
listed in Table 3.1 and illustrated in Figure 3.1 [23]. The effects 

of leakage through these and other data windows are illustrated by 


Ramirez [76]. 
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3.2 The Discrete Power Spectrum 

In practice, spectral estimates calculated according to the 
above procedure are produced with the aid of a digital computer. 
Therefore, a continuous signal must be expressed as a series of discrete 
values which generally represent samples of the signal amplitude 
collected at evenly spaced time intervals. If the process x(t) is 
sampled in such a manner and Fourier transformed as in Equation (3.2), 


the result is 


wipe ls ees n 
X (f£) = i 7 X(£-7 (3.8) 
n=—© 
where A is the time (in seconds) between samples [75]. If X(f£) = O when 
ere a, then l 1 
A x) Osea, 
X (£) = (GigS)) 
i Logegady sheseee 
A A Di ok 


This means that the Fourier transform of a process x(t) may be 
separated from the Fourier transform of the sample values whenever the 
magnitude of X(f) is zero for all frequencies greater than the so-called 
Nyquist frequency which is equal to 1/2A Hz. 

If the signal is not band-limited in this way, components 
with frequencies higher than the Nyquist frequency will appear at some 
frequency below the Nyquist frequency because of the summation in 


Equation (3.8). This effect is referred to as aliasing. 


3.2.1 The Discrete Fourier Transform 
If X (£) is to be evaluated digitally, Equation (3.2) cannot 
= 
be applied directly. Instead, the discrete Fourier transform is 


employed. The discrete Fourier transform of a series of numbers 
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(3280) 


where k = 0.1 .4.,N=-1 and j =4f-1. The result of this operation is a 
series of N evenly spaced samples of X(f£) spanning the Sreenererance 
between 0 and 1/A Hz. Only the first N/2 + 1 points of this series are 
unique. As Equation (3.9) shows, the remaining points are the complex 
conjugates of the first half of the series. 

Signals encountered in practical situations are not ideally 
band-limited as the sampling theorem requires. If, however, the 
magnitude of the high frequency components decreases with increasing 
frequency, the effect of aliasing can be minimized by sampling the 
signal at a rate high enough to insure that the magnitudes of the com- 
ponents above the Nyquist frequency are negligible in comparison to the 
magnitudes of the low frequency components. This requirement can be 
satisfied for any signal by passing the signal through an electronic 
low-pass filter before it is sampled. The cut-off frequency of the low- 
pass filter and the sampling rate are selected to yield a Fourier 
transform which covers the entire frequency range of interest with a 
minimum of aliasing. 

If this procedure is not acceptable for a particular signal, 
a specific procedure which is applicable only to that signal may often 
be implemented in order to reduce aliasing. An example of such a 
procedure is the technique proposed by French and Holden [19] for 
sampling neuronal spike trains. The spike train, which is treated as a 
geries of delta functions, is digitally filtered in order to generate a 


series of alias-free, evenly-spaced sample values. 
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A further error in the evaluation of discrete Fourier 
transforms is introduced when the sample values of the continuous signal 
are digitized for computer processing. Each sample value must be 
represented by a discrete binary number (quantization level) with a 
value which matches that of the sample value as closely as possible. 

The error ape eae by this procedure, called quantization error, can 
be viewed as noise which has been added to the original signal. Bendat 
and Piersal [14] show that the RMS value of this quantization noise is 
equal to 0.29 times the minimum possible difference between two quan- 
tization levels. This means that the RMS signal-to-noise ratio is equal 
to m/0.29 where m is the number of possible quantization levels within 
the amplitude range covered by the signal. Therefore, the allowable 
Bera roeneiee ratio determines the number of quantization levels 
required for a particular signal. Tufts, Hersey and Mosier [78] have 
shown that the effects of quantization noise on the Fourier transform 
are negligible when 256 or more quantization levels are available. 

Evaluation of the discrete Fourier transform of a series of N 
points according to Equation (3.10) requires approximately n? arithmetic 
operations. If, however, a fast Fourier transform algorithm is used, 
the number of operations required can be reduced to approximately 
2N*log.N [79]. This reduction in computational effort is achieved by 
judiciously factoring Equation (3.10). Fast Fourier transform algo- 
rithms are applied to time series for which N is an integral power of 2, 
but algorithms for transforming time series of other lengths are avail- 
able if required [80]. Otherwise, a time series may be extended to the 


desired length by attaching a train of zeroes to either end [80]. 
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3,2.2 The Pertodogram 

The square of the magnitude of a discrete Fourier transform 
is called a periodogram, and when this quantity is scaled correctly, it 
provides an estimate of the power spectrum. If X(f,) is the discrete 
Fourier transform based on a series of N samples of the process x(t) 
which has been modified by a data window w(t) and extended to length N' 
with a train of zeroes, it can be shown [80] that the resulting 


modified periodogram is 


A N y) 
oe Gos sara Ey, Grit) 
where 
N'-1 
U = a %  w(qA) Ain sys Coal) 
q=0 Cte a eee Ne 


The expected value of the modified periodogram is 
approximately Soe (FR) and its variance depends upon the shape of the 


data window [81]. For the rectangular window, 


Aw 


var[$._(£,)] = 2G ) (3.13) 


k 
Windowing with a tapered data window tends to reduce the effective 
length of a data segment because equally important data points are 
weighted unequally, and as a result, the variance of the periodogram is 
slightly inflated. However, Durrani and Nightingale [82] evaluated the 
variances of periodograms resulting from the application of various 
reasonable data windows to Gaussian data, and they found that, even for 
a severely tapered data window, the variance is inflated by only about 


10% over the value predicted by Equation (3.13). 
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The expected value and variance of a periodogram calculated 
from a particular time series are not affected if the time series is 
extended with zeroes before the periodogram is evaluated [80]. However, 
the correlation between adjacent periodogram ordinates increases as more 
zeroes are added to the time series [80]. The addition of extra zeroes 
reduces the spacing between periodogram ordinates without changing the 
shape of the spectral window. Jenkins and Watts [23] show that the 
covariance between two spectral estimates at different: frequencies is 
proportional to the amount of overlap of the spectral windows centred at 
these frequencies. The term "spectral window'' is normally associated 
with spectral estimates derived from the Fourier transform of the auto- 


covariance. function, For this situation 


2 
a pi a (3.14) 


where W, (f£) is the spectral window and Wy) is the frequency window 


Pears 


3.2.3 Smoothed Spectral Estimates 
A smoothed spectral estimate may be evaluated by averaging 
several modified periodograms (Equation (3.11)) from the same 


stationary process [84] 


na it 
Se SE) Sia So ey Cees) 
i=1 i 
The expected value of the smoothed spectral estimate [84] is 
x 2 
E{S_(£,)] =f H(f) S_ (£-£,)dé (3.16) 
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where See FD is the true power spectrum of x(t) and 
| N'-1 . 
i 
H(f) me | ot w(ahjal ea (3517) 
N'U 
q=0 
Thus the spectral estimate at a particular frequency is influenced 
(biased) by the power spectrum at other frequencies. The effect of this 
bias may be viewed as a decrease in spectral resolution, and its extent 


is dependent upon the width of the frequency band influenced by H(f). 


The variance of this smoothed spectral estimate [84] is 


| 


Var[S.(£,)1= ; Var[S,_(f,)]. (3.18) 


Therefore, if enough data is available, the variance can be reduced as 
much as desired simply by averaging over the required number of 
periodograms, The evaluation of a spectral estimate from a fixed 
amount of data involves a trade-off between variance and spectral 
resolution, If many periodograms are calculated from short, non- 
overlapping data segments and averaged, the variance of the resulting 
spectral estimate will be low. However, because the periodograms are 
calculated from short data segments, only a few sample values of the 
power spectrum will result, and therefore they will be widely spaced 
over the frequency band between 0 Hz and the Nyquist frequency. On the 
other hand, if periodograms are calculated from longer data segments in 
order to improve the spectral resolution, fewer periodograms will be 
available and the variance of the estimate will be higher. 

Welch [84] has shown that a decrease in the variance of a 
spectral estimate evaluated from a fixed amount of data may be achieved 
by averaging periodograms calculated from overlapping data segments. 


Such overlapping will increase the number of periodograms available 
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from a given amount of data, but, because the segments are not 
independent, the reduction in variance will not be as great as the 
reduction which could be achieved by averaging periodograms calculated 
from the same number of independent segments. The actual variance 
reduction for a particular overlap depends upon the shape of the data 
window, but Welch has suggested that if the segments overlap by half 
their length, a near maximum reduction in variance can be achieved for 
any particular window. For a data window similar to the Hanning window 
and an overlap of 50% of the segment length, the variance of the spect- 
ral estimate is reduced by approximately 40%. If periodograms from the 
same number of independent segments were averaged, the reduction in 
variance would be 50%. 

The multiplication of a time series by a data window or the 
convolution of its Fourier transform with a frequency window is referred , 
to as linear modification. A similar procedure, called quadratic mod- 
ification, may be applied to the autocovariance function or to its 
Fourier transform, the periodogram. With quadratic modification, a 
smoothed spectral estimate may be calculated by convolving a periodo- 
gram with a function known as a spectral window [23]. This is equiv- 
alent to multiplying the autocovariance function by a so-called lag 
window, the inverse Fourier transform of the spectral window. When 
spectral estimates are calculated from the autocovariance function, 
smoothing is accomplished by weighting the autocovariance function with 
lag windows similar to the functions used as data windows [27]. The 
resulting spectral windows are similar to the frequency windows 


illustrated in Figure 3.1. Spectral windows with much simpler shapes, 
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generally triangular, are now most often used when smoothing is carried 
out by means of convolution in the frequency domain [81]. 

The expected value of a spectral ree smoothed by 

quadratic modification is [23] 

1/2 
| = ss Ee ee tes Le) 
where W, (£) is the spectral window. This expression is.equivalent to 
Equatd.on. (34.16)» af W, (£) = H(f). 

The variance of a quadratically modified spectral estimate 

depends upon the shape and width of the lag window, w, (t) p23 
a Re cr 
var[(S (£1 © a Sy oae (3.20) 
where T is the total record length, 

Therefore, if the variance is to be low, the lag window must 
be narrow, but a narrow lag window corresponds to a wide spectral 
window. Since multiplication of the autocovariance function by a lag 
window is equivalent to convolution of the power spectrum by the cor- 
responding spectral window, this means that for a fixed amount of data, 
a trade-off exists between variance and resolution similar to the trade- 
off encountered in spectral estimates smoothed by averaging. 

Sloan [83], who compared the effects of linear and quadratic 
modification on a single periodogram derived from Gaussian data, found 
that linear modification results in a decrease in spectral resolution 
without a compensating decrease in variance, whereas the decrease in 
resolution due to quadratic modification is compensated by a comparable 


decrease in variance. He does concede, however, that linear 


cu 


ieaale ae : 7 


pees soe Sent 
nae Oey snag | 
| | ies 7 ) 
ies iii 
oo sbticsins et NOLAe2TgRS ant ‘wot too ata at ( 819 7 i 


nd Layne on ; eae on dof nti i 
sdamtiea lersseqe bolthbom vi Heol seybnup 8 a 7 sonatiey sar : 
{ec} (3) a Wobniw set vids lo dibiw ste aay 7a} “ogi 


a 
¢os.t) hits) a a he 0) jelasv 


®t fal 
edtansl liao94 (S203 offi al ¢ sisdw 
7 
; —-. a. i 
jean wobolw get ont ;wol ad-o% al acnstyay Sie IE (sot tert * 
Jaxs5eqe sbiw Bd ebnoqesti65 wohntw: gut wowed B ded swore od 


gelie vd qobionpl soo8 ts eyenasie ont he notisstiqhstund sorte poe 


@ 
-yo> sift vd murd59q2t9wng ‘sn4 Jo alban voenene ae tagtavinps eb w 


ysie® da. 3auoms boxtl & 10? sang Snmea elils evo sys 99qe acmml 
-ohe7) ewld oF taliwefe a it subdacry brit weeabawe To atone . 


‘duly eraves va bariapeme astknliate Co tba) nt borsi3myesns 
shiethanp bre sesail to eben aia rie ee Hy 
buypd Web, na bacuaid fond oe oa 6 MO 8 
botjitoass Petisaqe nb, a i estueor ne 169! 


D 
Hh 


modification is useful in reducing leakage when applied to data which is 
not Gaussian; and Otnes and Enochson [85] state that, for many practical 
situations, linear modification is a useful procedure because of its 
leakage-suppression properties. 

Jenkins and Watts [23] show that, for Gaussian white-noise 
processes, the probability distribution of a smoothed spectral estimate 
is proportional to a chi-squared distribution. The number of degrees 
of freedom for the distribution is 
2{E[S, (£,) 1} 


ae ie (327-45) 
var{S,_(£,,) } 


Vv = 


Thus the number of degrees of freedom of the approximating chi-squared 
distribution depends upon the variance of the spectral estimate which 
in turn depends.upon the smoothing procedure. Olshen [86] has shown 
that such a chi-squared approximation may be applied to a spectral 
estimate from any process representable as a moving average of 
independent (finite variance) random variables. Based on this chi- 
squared distribution, an approximate confidence interval for a 


particular spectral estimate may be calculated [23]. 


3.3 The Cross Spectrum 

The relationship between components of the same frequency in 
two signals determines the value of the cross-spectrum at that frequency 
[23]. If x(t) and y(t) are two stationary processes, their cross 


covariance function, ee is 


T/2 
Jim. f 


Bee (0) 2 rns op ip x(t)y(t+t) dt (3.22) 
xy =)? 
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The cross spectrum is the Fourier transform of R(T) 
xy 


-j27£T 


Bagh ed ey. Re dt | (323) 


Procedures for estimating the cross spectrum between two signals are 
similar to those for estimating the power spectrum of a single signal 
[81]. Fourier transforms X, (£) and Yp (fy) are evaluated according to 
Equation (3.10) from time series of length N which represent trains of 
simultaneous samples of the two signals x(t) and y(t). The cross 
periodogram Sree is Aree by 

Sale G 200 Mek galb si @ seo i 4 sh9 8 Nea (3.24) 

See oc U ee ok Tete 
where Xp (£,) is the complex conjugate of X(f,) and U is given by 
Equation (3.12). A smooth cross spectrum estimate may be obtained by 
averaging several such periodograms or by convolving a single cross 
periodogram with a suitable spectral window. The expected value of a 
cross spectrum estimate is approximately ye and the accuracy of 
this approximation improves if the width of the data window is 
increased [23]. 

In contrast to the power spectrum, which is always a real 
quantity, the cross spectrum at a particular frequency must generally 
be expressed as a complex number 
ge) - My (Eye eay Ae? (3.25) 


“a 


where ay) is the magnitude of the complex number and Sets) is its 


phase angle. For unmodified smoothed cross spectrum estimates [23774 | 


R il! 2} ald 
Van[M ei) = on es) : + seal (nO) 
Vey k 
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and 


var[Q,. (£,)] ® alles (3.27) 


where n is the number of cross periodograms averaged to produce the 
smoothed cross spectrum estimated, ey is the true cross spectrum 


: OE 
magnitude and es is a quantity known as coherence given by [14] 


pe gre Ey SOUTES (3,28) 


The number of periodograms which must be averaged in order to achieve a 
desired variance therefore, depends to a large extent upon the value of 


the coherence, a quantity which always falls between O and 1. 


3.4 Coherence 
The linear correlation between two processes x(t) and y(t) at 
: ; 2 
frequency f is given by the coherence, ee A coherence of 1 occurs 
for situations where x(t) and y(t) are linearly related, and a coherence 
of zero indicates that the two processes are completely unrelated. A 
coherence greater than zero but less than unity indicates that the 
apparent relationship between the two processes is nonlinear either 
because a truly nonlinear relationship exists or because the relation- 
ship is contaminated by extraneous noise [14]. The coherence value 
provides a measure of the degree of this apparent nonlinearity. 
Because the value of a coherence estimate based on unsmoothed 
spectral estimates will always be unity no matter what the true coherence 
: R2. 
is [9], a smoothed coherence estimate, Yyties must be evaluated 


according to 
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Benignus [87] experimentally evaluated confidence intervals for such a 
coherence estimate, and his results were confirmed by Carter et al. [88]. 
Benignus [89] also demonstrated that the confidence intervals are not 
greatly dependent on the probability distribution for the population of 
sample Seer in calculating the coherence estimate. Carter et al. 
[88] showed that the bias for coherence estimates increases as the true 
coherence approaches zero and as the number of degrees of freedom 
associated with the spectral estimates decreases. When n, the number 


of periodograms averaged, is large (say greater than 30) the bias is 
R2 1 2 2 
Be a eee (720) 


The variance of the estimate under similar conditions is 


i 2 ve 
var [y2 (£ )] s ne ete 0 
ae 2 Gab 
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Plots showing the exact Page and variance for various values of n and 
Ve tt) are presented by Carter et al. [88]. 

Theoretically, the coherence between two signals is not 
affected if the signals are subjected to linear filtering operations 
[90]. As Foster and Guinzy [90] demonstrate, however, this statement 
only applies to the sample coherence neeinea in Equation (3.29) if the 
filter transfer functions are identical or slowly varying over the 
length of the spectral window. If u(t) and v(t) are filtered versions 


of x(t) and y(t) where the respective filter transfer functions are 


A(f£) and B(f), then 
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W, (£) is the spectral window, and A*(f) denotes the complex conjugate 


Of AE); 


3.5 The Transfer Function Estimate 

If the two processes x(t) and y(t) are related by a linear 
system with a transfer function H(f), the following relationships apply 
[9,91]: 


ee = H(f) +S (£) (34,33) 


2 
Sy Cf) = |H(£)| “8. (£) (3.34) 


oN 


Based on these relationships, H(f,), an estimate for the transfer 


function, may be calculated from estimates of the auto and cross 


spectra. a 
r EGE y 
H(to)) =. G35) 
5 ena) 
xx k 
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H(f,) is a complex number which may be expressed as 


Niaeig tea APE) 
H(f,) = G(E, Je k (3.36) 
where 
F M (£,) 
Meg yer (370) 
fe een) 
xx k 
and 
P(f,) = ae (3.38) 


G(£) represents the gain of the transfer function, and P(f) describes 
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the phase spectrum. Approximate confidence limits for gain and phase 
are given by Jenkins and Watts [23]. | 

If the relationship between x(t) and y(t) is nonlinear, the 
quantity H(f) is the linear operator that best approximates y(t) in a 
least squares sense by acting on x(t) [13]. The coherence for such a 
system can be interpreted as a measure of the degree to which each 
frequency component of y(t) is related to the corresponding component 


of x(t) by a linear, time-invariant system, 


36 


ag SU ets i ay es iy: ey 
7 an ea 5 (1)¥ bows ix ‘geswitod geiplisictated aa t 
8 a fox aa demtsnaqg dapd tens norayego teoat! 
e doce tot snnaxsilod aft .ff1) C4) no ‘gntios vd senoe 


done dokdw o3 aevgph od3 Jo saueeom 6 28, bossiqradok ad nso a 


Inemodiies gnibneqesti02 sii. oF betelay at or ta sa vonsut 


.eedgeve doslbraval=s ; s 
aban reoakt wyd ¢ ad 


CHAPTER 4 
SPECTRAL ANALYSIS APPLIED TO SINGLE MOTOR UNITS 
4.1 Input-Output Relationship 

One approach to the development of a mathematical model for a 
physical system requires a set of experimental measurements of the 
system's input-output characteristics. In Many cases, especially for 
complex systems, such measurements provide the only available basis for 
a model. It is therefore important that the experimental procedure be 
designed to reveal the response of the system under conditions 
which are as natural as possible. 

The axon and muscle fibres of a single motor unit are com- 
ponents of a system which produces tension in response to electrical 
stimulation. In the body, the input signals to this system are trains 
of narrow spikes originating at the cell body of the motoneuron. 
Experimental measurements which are to be used in modeling the system 
should therefore be obtained for input signals of this form. Although it 
may be possible to develop a model for the response to signals other 
than spike trains, such a model would be of little use in illustrating 
the physiological behavior of the system. 

The tension developed by a group of muscle fibres depends 
upon such factors as stimulus rate [54,72], the length of the muscle 
[54,93,94] and the rate at which the length is changing [95]. The 
tension can also vary with time because of potentiation effects [72] or 
muscle fatigue [96,97]. Measurements from an experiment in which all of 
these factors are simultaneously varied would be difficult to interpret. 
A more reasonable approach would be to study the effects of varying one 


factor at a time. Measurements of isometric tension are particularly 
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suited to such an approach because the muscle length is fixed. Since 
the effects of fatigue or potentiation may be minimized by keeping the 
period of stimulation short, it should be possible to develop a model 
illustrating the effects on isometric muscle tension of stimulating the 
axon with a train of narrow pulses, 

A sigmoid-shaped curve relating isometric muscle tension to 
stimulus rate results when a group of motor nerve fibres [54,93] is 
stimulated with a series of evenly-spaced pulses. Although the slope of 
this curve is higher for intermediate stimulus rates than for low or 
high rates, the response to small fluctuations in stimulus rate is 
approximately linear. Presumably, a similar linear relationship applies 
for single motor units. If so, the system's frequency response may be 
determined for a variety of muscle lengths and mean stimulus rates, and 
a set of linear models relating changes in stimulus rate to changes in 
muscle tension may be developed. 

Muscles respond well to signals which have components at the 
low frequencies associated with normal movement [12]. However, the 
response to components with frequencies above about 8 Hz is weak, and 
for higher frequencies, the response is weaker still. If single motor 
units respond similarly, their overall behavior would be well- 
represented by frequency response plots covering a band of frequencies 
with some upper limit. 

An approximate value for this limit can be determined from a 
study of the motor unit's twitch response. A nerve action potential and 
the resulting isometric muscle tension for a typical motor unit are 
illustrated in Figure 4.1. The magnitude of the transfer function (gain) 


calculated from the Fourier transforms of these two waveforms is 
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FIGURE 4.1 Nerve Action Potential and Isometric Twitch Waveforms 
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illustrated in Figure 4.2. From this plot, which has been normalized 

so that the gain approaches unity at low frequencies, it can be seen 
that the gain drops off rapidly as frequency increases, For frequencies 
greater than 20 Hz, the gain is less than 3% of the low frequency gain. 
Therefore, 20 Hz appears to be a reasonable upper limit for modeling 


the frequency response of a single motor unit, 


4.2 Input Signals 

The conventional approach to the measurement of a system 
transfer function involves recording the responses to a series of 
Perecantroeduenor signals applied at the input. However, this approach 
cannot be applied to the study of single motor units because a pulse 
train containing only a single frequency component cannot be created. 
Therefore, a means for separating the Peeeueney components of both the 
input pulse train and the output tension record is required. The 
spectral analysis techniques described in Chapter 3 can be used for 
this purpose. If the system is linear for the input signals of 
interest, the resulting transfer function can be used to predict the 
response to any pulse train. If the system is nonlinear or noisy, the 
coherence estimate will indicate the degree of these effects and the 
transfer function estimate will provide the best linear approximation 
for the output in response to a train of pulses. 

Blackman and Tukey [73] state that the power spectrum of a 
signal which is to undergo spectral analysis should be reasonably flat. 
This criterion is quite well satisfied so long as the power spectrum 
does not change appreciably over the width of the spectral window. 


This means that the power spectrum of the input signal should not 
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contain any sharp peaks or valleys. A train of pulses with interpulse 
intervals which are periodically related is therefore not a desirable 
input signal. Bayly [100] has demonstrated that the power spectrum of 
such a pulse train would contain sharp peaks at frequencies dependent 
upon the periodically related interpulse intervals. A pulse train with 
random interpulse intervals is therefore more suitable as an input 
signal for the measurement of the motor unit transfer function. 


However, the probability distribution of the interpulse intervals must 
be chosen so that a pulse train irk a reasonably flat power spectrum 
results. 

The class of distributions known as gamma distributions can 


be used in this situation. For such a distribution 


g&,.8-1,-Bt 
'(g) 


P(t) = (4.1) 


where p(t) is the probability of an interpulse interval of time t, g is 
the order of the distribution, 8 is a scale factor, and ['(g) is the 
gamma function of g. 

Nelsen [98] has shown that the power spectrum of a train of 
impulses with interpulse intervals which are independently distributed 


is 
2 
1 Yalu (ft) | 
Se Cf) Salida eed ee a (4.2) 
|1—mM (£) | 
Ct 
where t is the average interpulse interval and M, (£) is the 


characteristic function corresponding to the probability distribution 
p(t) 
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At f=0 the power spectrum includes an impulse of size Wea If the 

0 Hz component is subtracted before the signal is Fourier transformed, 
this impulse will have no effect on the results of the spectral 
analysis. 

The plots in Figure 4.3 of the power spectra evaluated 
according to Equation (4.3) for various gamma distributions illustrate 
that for low values of g, the spectra vary slowly with frequency. In 
fact, when g=l1, the power spectrum is flat. A first-order gamma 
distribution is equivalent to an exponential distribution, and the power 
spectrum for such a distribution has been evaluated by Rice [99] for 
the case where the pulses have finite width. Neglecting the 9 Hz 


component, this power spectrum is 
eo) eC!) (4.4) 


where Se is the power spectrum of a single pulse and 8 is the 
scale factor for the gamma distribution. 

The input signal applied to the motoneuron is transmitted to 
the muscle fibres as a series of action potentials. Therefore, the 
power spectrum of each input pulse should be similar to that of an 
action potential. A train of narrow, rectangular pulses can be 
produced with little difficulty, and if the width of these pulses is 
chosen correctly, the difference between the power spectrum of the 
pulse and that of an action potential is negligible at low frequencies. 
The power spectrum of a typical nerve action potential is flat within 
1% over the frequency range 0-20 Hz, and the power spectrum of a 


3 msec rectangular pulse is a (sin x)/x function which decreases 


monotonically to about 99% of its low frequency value at a frequency 
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of 20 Hz. Therefore, if the transfer function is to be estimated over 
a frequency range of 0-20 Hz, a train of 3 msec rectangular pulses is a 
suitable test signal, because the difference between the two power 
spectra over this frequency range is never greater than 1%. 

A frequency response plot based on the use of such an input 
signal must be scaled according to the pulse amplitude. The scaling is 
required because the magnitude of a transmitted action potential is 
independent of the size of the pulse which initiated it; and, since the 
gain is measured with respect to the input pulse train, the measured 
gain will be inversely proportional to the magnitude of the stimulating 


pulses. 


4.3 Spectral Analysis Procedure 

The development of a specific procedure for the measurement 
of motor unit transfer functions by means of spectral analysis must 
include a certain amount of trial and error adjustment. An initial 
procedure, based on the available knowledge of the motor unit's freq- 
uency response may be devised, but the results of applying this proced- 
ure will probably indicate the need for modifications. Any such 
procedure must, of course, work within the limits of the available 
instrumentation and computing facilities; and, if possible, the results 


of each transfer function measurement should be immediately available. 


Because of its greater speed, the direct method for spectral 
estimation is preferable to the indirect method in this application. 
With the direct method, smoothing can be achieved either by averaging 
several consecutive periodograms or by convolving a single long per- 


iodogram with a suitable spectral window. Since the calculation of a 
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long periodogram with the fast Fourier transform requires more computer 
core memory than does the calculation of several short periodogranms, 
the averaging approach is preferable. If a means of displaying the 
spectral estimate is available, this display may be continually updated 
as more periodograms are added to the average, and the decision to stop 
sampling may then be based on the smoothness of the displayed estimate. 

The maximum frequency for which the transfer function 
estimate is to be evaluated determines the required sampling rate. 
Since the desired input signal contains significant components at 
frequencies above the desired Nyquist frequency, some form of low-pass 
filtering must be used to prevent aliasing. 

The method of French and Holden [19] would be suitable for 
eliminating this aliasing in the input pulse train, but this method 
requires special hardware to accurately measure the time of occurrence 
of each pulse. When this hardware is not available, an electronic low- 
pass filter may be used. Koles[9] has shown that a filter consisting 
of two cascaded second-order filters with Butterworth characteristics 
is eed Pa eater this application. If the corner frequencies of the 
Butterworth filters are equal to 60% of the Nyquist frequency, the 
aliasing error for a signal with a flat power spectrum is no higher 
than 6% at all frequencies below 80% of the Nyquist frequency. There- 
fore, a Nyquist frequency of 25 Hz and a filter corner frequency of 
15 Hz are suitable for measuring the frequency response of a motor 
unit for frequencies below 20 Hz. 

Only a limited amount of time is available for sampling the 
response of a motor unit subjected to continuous stimulation. Since 


motor units fatigue when thay are repetitively stimulated, the total 
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sampling time should be shorter than the time required for the effects 
of fatigue to become apparent; but, from the results of Smith and 
Lannergren [72], it appears that this time can be as short as 3 seconds. 
With a Nyquist frequency of 25 Hz, only 150 samples can be collected 
during this period of time, and if several periodograms are to be 
calculated from these samples, the frequency resolution will be very 
low. However, because the motor unit quickly recovers from the effects 
of a short burst of Semi one the frequency resolution may be 
improved by averaging individual periodograms calculated from the 
response to widely separated short bursts of stimulation. With a 
stimulus burst 2.56 seconds in length and a sampling interval of 
20 msec, independent periodogram coefficients may be calculated for 65 
frequencies spanning the range 0-25 Hz. 

When the power spectrum of a signal is reasonably flat over 
the width of the spectral window, the effects of leakage will be simi- 


lar for all periodogram coefficients. The relative magnitudes of the 


coefficients will therefore be unaffected. If an input signal with a 
flat power spectrum is used in studying the responses of single motor 
units, the effects of leakage between nearby coefficients will not be 
apparent. However, because the system severely attenuates the high 
frequency components of the input signal, leakage between widely 
spaced periodogram coefficients could noticeably affect the spectral 
estimate for the output signal. The Hanning window, with its wide 
major lobe and very small minor lobes therefore appears to be a suit- 
able choice as a data window. However, if the Fourier transformed 
samples are to be stored for further analysis, the rectangular data 


window is the most appropriate choice. With this approach, the stored 
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values may be subsequently convolved with any desired frequency window 


to reduce leakage. 


4.4 Summary 


The above discussion may be summarized as follows: 
The portion of a motor unit consisting of the axon coupled to the 
muscle fibres may be regarded as a single-input, single-output 
system. 
Although this system could conceivably respond to a variety of 
input signals, any model of the system's tm vtvo behavior should 
be concerned only with the response to trains of narrow pulses. 
The relationship between stimulus rate and muscle tension is 
influenced by many factors, and this relationship is, in general, 
nonlinear. However, a set of linear models fotecine isometric 
tension to the low frequency components of the input spike train 
may possibly be developed. 
Such models would be derived from plots of the system's low freq- 
uency response. Spectral analysis techniques may be used for such 
frequency response measurements. 
A suitable input signal for such studies is a train of narrow, 
rectangular pulses. The probability distribution of the interpulse 


intervals can be similar to that of a low-order gamma distribution. 
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CHAPTER 5 
PROCEDURE 

5.1 The Nerve-Muscle Preparation 

All experiments were performed on iliofibularis muscles taken 
from adult female specimens of Xenopus laevta. One end of this long, 
thin muscle which runs along the back of the thigh, is attached to the 
pelvis at the ilium, and the other end is joined to both the lower end 
of the femur and the upper end of the tibio-fibula by means of a split 
tendon. The muscle contains approximately 600 fibres of various types 
arranged into distinct zones. The lengths of the muscles used in the 
experiments, measured with the hind limb in a normal, relaxed position, 
ranged from 28-35 mm. 

The portion of the animal including the pelvis and the hind 
limbs was separated from the rest of the body, skinned and placed in a 
dissecting tray containing enough Ringer's solution* to cover it. The 
iliofibularis muscle in one limb was exposed, and its length was 
measured with the limb in a natural relaxed position. The muscle was 
then removed from the limb together with the portion of sciatic nerve 
approximately 4 cm in length beginning at the sciatic plexus and running 
to a point approximately 1 mm distal from the exit of the nerve branch 
serving the iliofibularis. This dissection was performed carefully to 
insure that the nerve branch entering the muscle was not stretched or 
otherwise damaged. The nerve-muscle preparation was then placed in a 


shallow plexiglass container filled with Ringer's solution. A small 


*The Ringer's solution had the following composition (mM): NaC1-112.0, 
KCI-2.5, CaCl9-2.0, NaHC03-2.5, and glucose-5.5. The solution was 
saturated with a mixture of 95% O09 and 5% COo. 
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ete was made in the tendon at each end of the muscle, and the muscle 
was fixed between two stainless steel hooks inserted through these holes, 
The nerve branch was supported on a thin plexiglass ledge fastened to 
the bottom of the container, The connective tissue sheath covering the 
sciatic nerve was slit longitudinally at the distal end over a distance 
of about 3 mm and peeled back in order to expose the cut ends of the 
nerve fibres within. All of these ends were then carefully pulled back 
until only the intact fibres which entered the iliofibularis remained, 
and a dissecting needle was used to cut all but one of these intact 
fibres. 

A small plexiglass block with a concave depression in the top 
was placed in the container with the nerve-muscle preparation. The 
proximal end of the sciatic nerve was placed in the depression, and the 
depression was filled with Ringer's solution. A seal of petroleum jelly 
was placed across the nerve and the plexiglass block immediately distal 
to the point where the nerve entered the depression in the block. The 
whole container was fastened to a small table inside a Faraday cage, and 
the wire hook passing through one end of the muscle was replaced by a 
thin stainless steel wire hook attached to a Devices Instruments type 
2STO2 isometric force transducer. The resonant frequency of the force 
transducer and the attached wire hook was approximately 250 Hz, and the 
compliance of this arrangement was about 0.01 mm/g-force., A silver- 
silver chloride electrode was lowered into the pool of Ringer's 
solution in the container, and a similar electrode was lowered into the 
pool of Ringer's solution contained within the depression in the plexi- 
glass block. The completed preparation, itlisttated in Figure 5.2. 


was at room temperature (21°-23°C) throughout the experiment. 


vesfor enact Mga’ borveant alaon fesse es ows ‘neewied 6 
07 bensteai sxbel aenlgiralt nidd a a0 betaoggue pape 
sda gnizevoo fteade gueati svistaennoo dT “erantna09 aria ‘to wodsed eda 
eorgsath 9 r8v0 bas !esath of) 38 yffentbustgnol alle sew oviad obtatoe 
soda 2%. has tug 913 seoqea oF q4br6 nk Yond belesg bae mm € 200de 20 a 
iopd bollug Ptutayas Abdi sirsw, aay seats ip LA -abitaiw aesdes ownage ; 
,beatemss etyaiudtiet!t 211 bereics dostdw aaxdl? tontak off, uno Lisa . 
tonppk sasns to eno tud ffs Jus od baam ve afbaon agiiosserb a bas . J 

qos sis mi satiawornah oveoncs a Miiw dscld sesiptesla ofeme A : 
sav .notvetsqesig s{oeum-avier sit atv tantsangs orl ink beoalq eBw 


ef? baw ,cobeesrqsd offi at beoslq eaw sviso sitetse ody To bas Lemkxorq 


yilet ausiortsg Io lose A .notjutoe e!q7agntt dalw belli? sew nobeeezgeb 
Ietath yloistbaant anotd eaaigtxstq ona “par “ovien sd9 s8otss beseiq aay 
efT .doold aft af noleesrqeh oft chai avian sdv svedw Integ dg 03 
ban ,9gao yebeanT o sblent glds2 bana: 8 ot beaases2 @aw tantssnes elodw : 
a y@ beoglass exw oloenm oft to bes Sno dauoxds antaneg, acon sxtw eds "4 | 
sqyd aigsmyrienl agntvsd a ot badossse sxond sale laste eesintase atdy 
sor0t adit to yoasuper? Jorsoesy oAT .rs9vbanats sod oksgeupat SoTeS 7 
ait bra ,sH O28 ylevemixoxrgqs. esw ool sitw barosdie adj bas Yoopbeasx3 di 
-1s7[tea A ent0t=—\n £f,0 juode asw Inemegmex2e ais Yo sonehfqmos ne ) 
e'tagnrh to loog ofd otek boitswol esw aborssafe aa _ 
aid ogni betevel exw shord26ls- Sarheee A bes Fanaa: ait: 
Pap pe ond mt saphignes iia niidiw bsnbetnos ENN, ee 


From Stimulator 


Plexiglass Chamber 


Plexiglass 
Block 


Support 


Single 
Nerve 
Fibre 


Isometric 
Force Transducer 


Petroleum Jelly 


See 
AZZ: 


Section A-~A 


SS 


Ses 


FIGURE 5.1 The Nerve-Muscle Preparation 
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The input signal was applied to the preparation by means of 
the two silver-silver chloride electrodes which were connected to one 
output of an Ortec 4710 dual-channel stimulator. Action potentials were 
generated in the motoneuron by applying a potential difference across 
the petroleum jelly seal with the stimulator, Such a potential 
difference caused current to flow along the axoplasm of the various 
nerve fibres and through the nerve cell membranes, thereby creating a 
potential difference across each membrane, When this potential 
difference reached the threshold potential, an action potential was 
generated and propagated along the axon, but only those action pot eri 


tials generated in the intact fibre could reach the muscle. 


Sekt Stimulation and Recording 

The stimulator consisted of two separate pulse generators 
triggered from a common source, The second output was connected 
through one of two identical low-pass filters described in Appendix 1 to 
one input of a 16-channel Hewlett-Packard 5610A analog-to-digital con- 
verter. If desired, the output of the stimulator could be gated with 
an externally applied signal. The stimulator could be triggered either 
by means of an internal clock or by an external signal. The external 
trigger signal was provided by a Precision Instruments PI-6200 tape 
recorder with an analog tape recording of a random pulse train. Several 
different pulse trains were recorded on the tape so that a pulse train 
with interpulse intervals distributed according to a particular 
probability distribution could be selected by positioning the tape at 


the desired point. The procedure for generating these pulse trains is 


described in Appendix 2. 
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The force transducer was connected to the input of a Tektronix 
3C66 carrier amplifier, and the output of this amplifier was passed 
through a low-pass filter to the second input of the analog-to-digital 
converter, The output of the carrier amplifier was also connected to 
one channel of a dual-channel Sanborn 320 chart recorder. The second 
channel of this recorder displayed the output of a frequency-to-voltage 
converter [101] which measured the instantaneous pulse rate at the 
stimulator output. An audio amplifier and loudspeaker were connected 
to the stimulator output so that this signal could be monitored aurally. 
The erranvenent of the equipment used in the experiments is illustrated 
in Figure 5.2 

Note that the input signal applied to the motor unit and the 
input signal sampled by the analog-to-digital converter originate from 
different outputs of the stimulator. This arrangement was adopted as a 
matter of convenience because the magnitude of the pulse required to 
initiate an action potential in the motoneuron varied with each 
preparation, The second channel of the stimulator produced a pulse at 
the same time as the first channel, but the magnitude of the pulse was 
not changed from experiment to experiment and, as a result, the 
magnitude-dependent scaling factor mentioned in Chapter 4 was the same 
for every experiment. In addition, the width of the stimulating pulse 
could be changed without affecting the results of the analysis. 

The threshold voltage for the preparation was considered to be 
the minimum stimulator voltage required to produce an observable con- 
traction of the motor unit muscle fibres. Those motor units which pro- 


duced a muscle twitch in response to a 1 msec stimulus pulse were 
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assumed to be fast motor units, Slow motor units were considered to be 
those which exhibited no twitch response but which responded with a 
slow muscular contraction to a train of 1 msec stimulus pulses at a rate 
of about 20 pulses per second. Any preparation which produced no 
observable single motor unit response was discarded. The intact nerve 
axon of such a npeparacics was assumed either to belong to a sensory 
neuron or to be damaged. 

Once the threshold voltage had been determined, the responses 
of the motor unit to a series of constant-rate pulse trains were 
recorded on the chart recorder, For these measurements and all suc- 
ceeding measurements, the stimulating pulse amplitude was set at 
approximately twice the measured threshold voltage. Generally, the 
stimulating voltage could be increased to a level at least 10 times 
higher than the threshold voltage before current spread resulted in the : 
initiation of action potentials in the cut ends of the axons leading to 
the muscle. 

Next, data required for the estimation of motor unit transfer 
functions were collected. The motor unit was subjected to periods of 
continuous stimulation, and its response to this stimulation was 
measured and recorded according to the procedure to be described below. 
Both the length of the stimulation periods and the specific stimulating 
signals employed varied from experiment to experiment, but, in most 
cases, periods of stimulation exceeding 1 min. in length were avoided. 
In an attempt to reduce the effects of muscle fatigue, the muscle was 
allowed to recover for several minutes after any prolonged period of 


stimulation. 
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5.3 Data Analysis During the Experiment 

Data collection for the evaluation of transfer function 
estimates was controlled by a Hewlett-Packard 2100S computer. The 
computer program which performed this data collection is illustrated by 
means of a flow chart in Appendix 3. The sampling times of the analog- 
to-digital converter could be controlled either by the computer or with 
a pulse train from a Hewlett-Packard 2781A pacer. In either case, a 
command to initiate sampling caused the analog-to-digital converter to 
sequentially sample the two signals applied to inputs 1 and 2 with a 
time interval between the two sequential samples of about 12 usec. 
Each sample was converted to a 10-bit binary number which was passed to 
the computer as part of a 16-bit binary number. The additional 6 bits 
carried control information. 

The sample values were collected as paired blocks of data. 
All of the points in one of these blocks represented samples from 
Channel 1 of the analog-to-digital converter, and all of the points in 
the other block represented samples from Channel 2. For each data 
block, the mean value of all the points was subtracted from every point 
in the block, the resulting series of values was Fourier transformed, 
and the results of this operation were stored on a disc memory. The 
collection of any data block pair and the Fourier transformation of the 
previously collected pair of data blocks proceeded simultaneously. In 
_ addition, aa average power spectrum estimate was calculated for the 
signal representing the motor unit output. This estimate was displayed 
on the oscilloscope as a plot of power versus frequency. The average 
power spectrum and the display were updated every time a new pair of 


data blocks were Fourier transformed. The smoothness of the spectral 
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estimate could be judged from the display, and the sampling procedure 
could be stopped whenever enough data blocks had been collected to 
yield an acceptably smooth estimate. Alternatively, a command to 
collect a fixed number of data blocks could be issued if desired. 

The length of the data blocks and the sampling rate were 
specified before sampling commenced. The number of points in any data 
block was limited to integral powers of two, and the maximum block size 
was 256 points. Samples could be collected continuously, or the pacer 
could be used in conjunction with the counter-controlled gate circuit 
described in Appendix 4to collect samples in bursts of a specified 
length. With this procedure, a large number of data blocks could be 
collected by subjecting the motor unit to many short periods of stim- 
ulation. When the motor unit was allowed to rest between periods of 
stimulation in this way, the effects of fatigue were reduced in com- 
parison to the case where all of the data blocks were collected during 
one long period of stimulation. 

The maximum possible sampling rate was limited because of the 
time required to calculate the Fourier transform. The collection time 
for a data block pair could not be any shorter than the time required to 
Fourier transform the previously collected pair of data blocks and write 
the results on the disc. With an efficient fixed-point fast Fourier 
transform subroutine [77] written in assembly language, data could be 
collected continuously on two channels at a maximum rate of 300 point 
pairs per second. 

When enough data blocks had been collected to yield adequat- 
ely smooth spectral estimates, transfer function estimates based on any 


portion of the recorded data could be calculated and displayed on the 
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oscilloscope as plots of gain, phase, or coherence versus frequency. A 
flow chart for the computer program which controlled the production of 
these plots is presented in Appendix 3. The values for gain in such 
plots indicated the gain as calculated from the filtered input and 
output signals. For all subsequent analysis, these values were mult- 
iplied by the gain setting of the strain gauge amplifier. The ae 
for gain expressed in this way are g-force/volt. Because the gain was 
inversely proportional to the magnitude of the input pulses, this 
magnitude was kept constant from experiment to experiment, as were the 
gains of the anti-aliasing filters. Therefore, any differences among 
the various gain plots were due only to differences in the responses of 
the motor, units. 

Because plots of the transfer function estimates were 
available immediately after each period of data collection, modificat- 
ions of the experimental procedures based on the results of each trans- 
fer function estimate could be carried out during the course of the 
experiment. After 12 preparations had been studied in this way, it 
was decided that, although such plots were useful in roughly gauging 
the responses of single motor units, further analysis, which could 
only be carried out after the completion of the experiment, would be 
necessary. Sampling intervals for these first experiments ranged 
from 10 to 100 msec, and mean stimulus rates ranged from 5 to 30 
stimuli/sec. In most cases, the muscle was kept at its natural 
resting length. The results of these experiments were used in 
determining a suitable experimental procedure for future studies. 

Samples for subsequent experiments were collected at a rate 


of 50 pairs per second in blocks of 128 sample pairs. Generally, 10 
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block pairs were collected for each transfer function estimate. The 
cut-off frequency for the anti-aliasing filters was 15 Hz. The 
Fourier transformed sample values for this group of experiments were 
stored on magnetic tape for further analysis by means of the computer 
program illustrated by a flow chart in Appendix 3. The results of 
212 spectral analysis studies for 34 preparations were stored in this 
way. 

Because the preparations were subject to fatigue, it was not 
possible to study each preparation with a full range of stimulus 
signals and muscle lengths. The effects of changes in muscle length 
were studied only for mean stimulus rates of 10 stimuli/sec. This 
rate was high enough so that the muscle tension waveform was not 
simply a series of individual twitches, and low enough so that fatigue 
effects did not appear rapidly. Various muscle lengths spanning a 
range of 2 mm above and below the natural resting length were studied 
in 17 of the 34 preparations referred to above. The effects of changes 
in mean stimulus rate were studied only with the muscle at its natural 
resting length. Stimulus signals having mean rates of 0.2, 5, 10, 15, 
and 20 stimuli/sec were used. Two or more of these signals were used 
on 22 preparations. Continuous stimulation was used on 14 preparations, 
and bursts of stimuli were used on 26 preparations. 

For 3 preparations, the time delay between the occurrence 
of the stimulus pulse and the beginning of the muscle twitch was 
measured by simultaneously displaying the twitch waveform and the 
stimulus pulse on a Tektronix 604 storage oscilloscope. When the 
displayed twitch responses were compared with the twitch response 


plotted in Figure 4.1, it was concluded that the beginning of the 
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twitch response could be determined to within about 0.5 msec from the 


oscilloscope display. 


5.4 Recording Action Currents in the Muscle 

In a few experiments, action currents associated with the 
contraction of the isolated motor unit were recorded. The level of 
the Ringer's solution covering the nerve-muscle preparation was 
lowered so that the surface of the muscle was exposed to the air. A 
layer of paraffin oil deep enough to cover the preparation was then 
floated onto the surface of the Ringer's solution, and two silver- 
silver chloride electrodes were placed in contact with the upper 
surface of the muscle. 

Action currents in the muscle produced a potential difference 
between the two electrodes. This potential difference was amplified 
with a Tektronix RM122 preamplifier and displayed on a Tektronix RM502A 
oscilloscope. A Grass Instruments Kymograph camera was used to record 


the oscilloscope display. 


5.5 Histological Techniques 

At the completion of each experiment, the external diameter 
of the intact axon was measured to the nearest 0.5 um with an ocular 
Micrometer . 

With several of the preparations, an attempt was made to 
identify the specific muscle fibres constituting the single motor unit. 
After all the transfer function measurements had been completed, the 
motor unit was fatigued by repeated stimulation. Alternate periods of 


rest and stimulation were continued until the motor unit produced no 
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observable tension when stimulated. This procedure was carried out with 
the expectation that it would cause depletion of the glycogen in the 
muscle fibres [103]. The muscle was then quickly frozen in a mixture 
of dry ice and methylbutane, and transverse sections approximately 8 um 
in thickness were cut with a freezing microtome. Serial sections were 
stained for succinic dehydrogenase [102] and for wiveeeen using the 
periodic acid-Schiff reaction [103]. In this way muscle fibres 
exhibiting glycogen depletion could be related to the muscle fibre 
types identified by Smith and Ovalle [34]. 

Another approach to relating motor unit types and muscle fibre 
types was attempted with several preparations. A single muscle fibre 
from the isolated motor unit was dissected out and prepared for electron 


microscopy according to the procedure used by Smith and Ovalle [34]. 


5.6 Mechanical Properties of the Whole Muscle 

A procedure: similar to that used by Halpern and Alpert [104] 
was used in studying the effects of the muscle's dynamic mechanical 
properties on the observed motor unit transfer function estimates. The 
iliofibularis muscle was removed from the animal and placed in the same 
plexiglass container used for the motor unit studies. One end of the 
muscle was connected to the force transducer and the other end was con- 
nected directly to the shaft of a Ling Altec Ltd. three ohm vibrator. 
The displacement of this vibrator shaft was sensed with a Hewlett- 
Packard 24DCDT-050 displacement transducer, and an electronic feedback 
circuit [9] insured that the gain of the entire puller system was flat 
within 1 db for frequencies below 100 Hz. A Gaussian white noise 


signal[9], band-limited to the range 0.4-80 Hz, was applied to the input 
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of the muscle puller, and transfer function estimates describing the 
relationship between the displacement of the puller shaft and the change 
in muscle tension were determined for the frequency range 0-100 Hz. 
These transfer function estimates were calculated with the same set of 


computer programs used in the motor unit studies. 
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CHAPTER 6 
ANALYTICAL METHODS 
6.1 Time Response 

In order to compare the isometric responses of the various 
motor unit types, several descriptive parameters were defined. The 
twitch amplitude was considered to be the maximum value of the muscle 
twitch tension resulting from a single stimulus, and the time between 
the initiation of tension increase and the occurrence of this peak 
tension was called the twitch peak time. Since the tetanic tensions 
for nearly all motor units did not increase when the stimulus rate was 
raised above 50 stimuli/sec, the maximum tetanic tension was defined 
as the peak tension produced in response to a burst of stimulation at 
a uniform rate of 50 stimuli/sec. In order to prevent fatigue damage 
to the motor unit while measuring this quantity, stimulation was 
terminated as soon as an obvious peak in tetanic tension had been 
reached. The ratio of twitch amplitude to maximum tetanic tension was 
called twitch/tetanus ratio. 

In studying the fatiguability of fast motor units, it was 
found that although all motor units could quickly recover from the 
effects of 10 sec of stimulation at a uniform rate of 30 stimuli/sec, 
such a period of stimulation was long enough to indicate differences 
in fatiguability. The fatigue coefficient was therefore defined as the 
tension at the end of this period of stimulation divided by the 
maximum tension produced during the 10 sec interval. 

Differences in the rate of tension development in slow motor 
units were most apparent at low stimulus rates, but tensions were low 


at such rates, and differences in tension were difficult to measure. 
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It was therefore decided that a rate of 20 stimuli/sec was the lowest 
suitable rate for studying tension development in slow motor units. 
The initial tension ratio was defined as the ratio of tension after 

1 sec to the tension after 20 sec of such stimulation. 

The tension and time values required to evaluate the above 
parameters were measured from the strip chart records collected during 
the experiments. The tensions recorded on the strip chart records 
represented the total muscle tension minus the resting tension. 
Therefore, all tension measurements indicated the change in isometric 


tension resulting from stimulation of the motor unit. 


6.2 Frequency Response 

Transfer function estimates were calculated from the data 
which had been recorded on magnetic tape during the experiments. These 
estimates were calculated by a Hewlett-Packard 2100 computer and 
plotted on a Tektronix 4010-1 display terminal. When possible, each 
measured transfer function estimate was then fitted with a curve rep- 
resenting the frequency response of a linear system. The computer 
program which controlled these operations is illustrated by means of 


a flow chart in Appendix 3. 


6.2. Choice, of Model. 

Initial examination of the plotted coherence estimates 
indicated that they could be divided into two categories related to 
motor unit type. For slow motor units, the points representing the 
coherence were scattered between 0 and 0.5 for all frequencies, 
thereby suggesting that the observed responses of these motor units 


would not be appropriately represented by linear models. For the 


a “4 A 
ada “steutove. a2) 
Nee nadaetion yer redo re veld mos) f 


‘= a | 
abroo9% vane une vl fi | barbed LAER r off i 
* peer day ; 


hat Lupe " ae 


+- 


idkensd goites7 oft uPrer olan ots Istox pa nf ; 
‘ 


ei 
jriteamoat ai sgnnis 3 soi oe eminent 
J 


sii x6d00 sda to io fists tie neti camelimiaitial’ 
, 7 
7 
sanoqeat eaaapeit a 


nteb sdj-mord toisluelan arswogatamivess fnka9ai sotto w os 


antvuh naa" oid enges 110. bobions oad bd ad Hi: 


> 


gasdt ceinsmiveqxs si 


bos wejponds OOLS hrndvatodaslbuee B ei fis batts id coe % las 
ane a 
dase ,efdleroqg ooriW tandiae yalqetb I-OL08 Ree no 


a, ‘ “on 
-~qaq7 svivo s ftiw beiil) inedy aew Stamh jes gubimcians sae? eT 
ous 3 : 


s93vgmos Sat .masaye saat 5 20 pentnea! oar ‘ona Eyl 


Yo ancom xd bossrteults et enokseveqe saont vettodiwes doiew nga 
2 xibasqqa wk sunty ¥ 


ra oa 


- 


pesumtses soneiadoy bys totq a Pons 
i 
od bess fer ae ah r= int) bebivib od tus D3 
vod sapanasotanyeatnicn ont deat ae & wo 
o ’ ~n i 
ee 163 ib bat red heres 
J 0 haw sat hi 


ars 


65 


fast motor units, the coherence was near some constant value (usually 
above 0.7) for all frequencies less than about 15 Hz, but above this 
frequency the coherence declined rapidly to near zero at 25 Hz. 

A study of the gain and phase plots for fast motor units 
revealed that the gain was constant for low frequencies and decreased 
at a rate of 12 db/octave at high frequencies. The phase approached 0° 
at low frequencies and was approximately -90° at the point where the 
low and high frequency asymptotes to the gain curve intersected. For 
high frequencies, the phase was in the neighborhood of -180°. The 
shape of these frequency response curves was therefore similar to the 
shape of the gain and phase curves for a second-order low-pass filter; 
SO, because of this similarity, such a system was chosen as a model 
for fast motor units. 

A time delay of about 10 msec between the initiation of a 
nerve stimulus and the beginning of the resulting muscle twitch could 
be measured directly, and the effects of such a delay were observed in 
the transfer function estimates as an increase in phase shift beyond 
the expected -180° at the high frequency end of the plot. Such a time 
delay was therefore combined with the second-order low-pass filter to 
produce frequency response curves which could be fitted to the measured 


transfer function estimates. 


6.2.2 Bias in Coherence Estimates 

An electronic low-pass filter with a response similar to that 
of a single motor unit was used to demonstrate the existence of bias in 
the coherence estimates. A coherence estimate for this electronic 


filter was calculated according to the same procedure used for the 
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Coherence 


1 10 


Frequency (Hz) 


FIGURE 6.1 Coherence Fstimates for Electronic Filter 


motor units, but the input signal was band-limited Gaussian white 
noise* [9] rather than a random pulse train. This estimate, illus- 
trated by the plotted points in Figure 6.1A, represents the coherence 
between a random signal passed through an anti-aliasing filter and the 


same signal passed through an identical anti-aliasing filter in series 


* The power spectrum of this signal was flat within 3 db from about 
0.04 Hz to 15 Hz. At high frequencies, the magnitude of the 
spectrum decreased at a rate of 20 db/octave. 
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with the electronic low-pass filter. The sharp drop in coherence at 
high frequencies is similar to the drop observed in the coherence 
estimates for single motor units when a rectangular data wandow was 
used, 

The effect of linear filtering on the coherence estimate 
can be determined with Equation (3.32). For ene situation described 
above, x(t)=y(t). Therefore Yay (£) should equal 1 at all frequencies. 
A(f) is the transfer function of the anti-aliasing filter, and B(f) is 
the transfer function of the electronic low-pass filter in series with 
the other anti-aliasing filter. The dotted line in Figure 6.1A shows 
the coherence estimate for this situation calculated with Equation 
(3.32) for a rectangular data window. The agreement between the 
dotted line and the plotted points indicates that the observed drop in 
coherence at high frequencies is largely due to bias caused by linear 
filtering [90]. 

Jenkins and Watts [23] state that bias in the coherence 
eeeee will increase with the length of the time delay between the 
two signals. The effect of the short time delay observed in single 
motor units was simulated by adding a 10 msec delay to the system 
described above, but the measured coherence estimate did not change 
noticeably, and the values for the theoretical estimate changed by 
less than 0.5%. These results indicate that the time delay associated 
with a single motor unit is too short to appreciably affect. its 
coherence estimate. 

The bias due to filtering can be reduced if leakage through 
the spectral window is decreased. This effect is illustrated in 


Figure 6.1B which shows the coherence estimate for the electronic 
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low-pass filter after the Fourier transforms of the input and output 


signals have been convolved with the Hanning frequency window. The 
dotted line shows that, in theory, the effects of bias due to filtering 
can be almost eliminated when the Hanning window is used. But, because 
of the greater width of the Hanning window's major lobe, leakage bet- 
ween adjacent frequencies is increased, Therefore, for this example, 
the coherence value nearest 0 Hz becomes slightly biased because of 
leabase from the 0 Hz component which has a value of zero. However, 
this slight bias at low frequencies is overshadowed by the improvement 
in the coherence estimate at high frequencies. The Hanning window is 
therefore highly suitable for reducing leakage in spectral estimates 
for single motor units. 

When the Hanning window is used, the coherence estimate 
still decreases slightly at high frequencies even though the theoret- 
ical bias is almost zero. However, two other possible sources of bias 
exist--aliased power and quantization noise. The gain of the anti- 
aliasing filter drops sharply at high frequencies, so the effects of 
aliasing will be most apparent near the Nyquist frequency. As well, 
because the electronic filter tends to reduce the magnitude of the 
high frequency components, quantization noise will be more apparent at 
high frequencies. Since band-limited white noise was used as the input 
signal to the electronic filter, the anti-aliasing filters could be 
removed without totally degrading the spectral estimates. If aliasing 
is the cause of bias in the coherence estimates, removal of the anti- 
aliasing filters should increase the bias. As Figure 6.1C illustrates, 
however, the bias is greatly decreased when the anti-aliasing filters 


are removed. At a frequency of 23 Hz, the aliased power in the input 
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signal is increased by a factor of about 10, and the aliased power in 
the low-pass filtered version of this signal is increased by a factor 
of about 6. Yet the bias at 23 Hz is about 10 times smaller without 
the anti-aliasing filters. However, the ratio of quantization noise to 
Signal amplitude at this frequency is decreased by a factor of about 10 
when the anti-aliasing filters are removed. 

The coherence between two linearly-related signals with power 
spectra X(f) and Y(f) is reduced when the signals are contaminated with 
added noise. If the spectra of the noise signals are N(£) and ey 
respectively, the coherence is [90] 


ep (6.1) 


a, N (£)77 N (£) 
| het 

When this expression is evaluated at 23 Hz for a signal which spans 
600 quantization levels at low frequency, the coherence is 0.956. In 
eos 6.1C, the coherence at 23 Hz for signals with low frequency 
amplitudes spanning about 600 quantization levels is 0.95. A close 
agreement between measured and theoretical ne is also apparent 
for other frequencies and for the situation where the anti-aliasing 
filters are included. It may therefore be concluded that the bias in 
the coherence estimate which remains after the signals are treated with 
the Hanning spectral window is due mainly to quantization error. 

Gain and phase estimates for the electronic low-pass filter 
determined with and without the Hanning spectral window are teeta 
in Figure 6.2. When these plots are compared with the plots of the 
corresponding coherence estimates in Figure 6.1, it becomes apparent 


that the gain and phase plots are noticeably affected only at 
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FIGURE 6.2 Transfer Function Estimates for Electronic Filter 
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frequencies for which the coherence estimate is severely biased. When 
the Hanning frequency window is used, the factors which cause the 
coherence estimate to become biased at high frequencies have little 


effect on the gain and phase estimates. 


6.2.3 The Second-Order Low-—Pass Filter 


The transfer function of a second-order low-pass filter is 


[105] 


k 


HUi£): = (6.2) 


l= payee = GREY AE 
n n 

where © is the damping ratio, and fe is the undamped natural frequency. 
The gain of the filter at a frequency of 0 Hz is equal to k. Since 
such filters are commonly used for electrical signals, k is often 
called the de (direct current) gain of the filter. Normalized plots of 
H(jf) versus frequency for several values of Zt are shown in Figure 6.3. 

For a given filter, approximate values for k, TZ, and f may 
be obtained from a plot of the filter's frequency response. The dc gain 
is indicated by the low frequency asymptote, and the natural frequency 
is defined as the frequency at which the high and low frequency 
asymptotes intersect, or as the frequency where the phase is equal to 
-90°. The damping ratio can be estimated by comparing the plot to a 
series of normalized frequency response plots such as those illustrated 
in Figure 6.3. Note that a low damping ratio results in a peak in the 
gain curve and a rapid change in phase near the natural frequency. 

The variability of such estimates for k, Tf, and fe ace 
upon the nature of the frequency response plot. This plot is based on 


measurements of the system's response at a number of different 
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frequencies. If the variability of these measured values is high, or 
if measurements are available for only a few frequencies, the positions 
of the best-fitting low and high frequency asymptotes become difficult 
to determine. Estimation of the damping ratio also becomes difficult 
under such conditions. Because of these oro: a more reliable 


method for fitting curves to transfer function estimates is needed. 


O24 Curve Fitting 
One scheme for fitting parameters to a transfer function 
estimate involves expressing the fitted transfer function, H(jf), as the 
ratio of two polynomials: 
py + Pp (if) + p,(4f)° +... 


' 7 (je) 
H(j£) =— = ; (6.3) 
laa ee Mapes a, (48)? on Ok, 


At a particular frequency f the difference between the transfer 


k? 


function estimate, H(f,); and the fitted transfer function, H(jf,), is 


e. where 


R P(jf,) 


e, = H(jf,) ~ UGE (6.4) 


k 


Levy [106] proposed a method of evaluating P(jf) and Q(jf) 
based on minimizing the sum over all experimental points of a weighted 
2 
squared error ley! where 


ey = e, Qf) (6.5) 


With this criterion, the coefficients of P(jf,) and Qcjf,) may be 
obtained as the solution to a set of simultaneous algebraic equations 


of the form 


[A][X] = [B] (6.6) 
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An approach similar to this was used with limited success by Kislack 


[107] for fitting parameters to transfer function estimates. 
A technique developed by Sanathanan and Koerner [108] allows 
the coefficients of P(jf,) and Q(if,) to be determined iteratively. 


With this procedure, the sum of en 2 is minimized where 


os Ole Qt, 


7 I ce ee 
k OCjE) S41 


EGan) 


and L is the iteration number. Since this expression converges 
rapidly to ene the weighting function associated with the error term is 
unity rather than Q¢jt,) as used by Levy [106]. Sanathanan and Koerner 
demonstrated that curves fitted on the basis of an unweighted error 
term follow the transfer function estimate of a complex system more 
closely than do curves fitted according to an error term which is 
weighted with Q(jf,). 

The iterative procedure has been applied by Streets and 
Hemphill [109] to the problem of fitting paraneceen to systems which 
include a time delay. With this approach, the time delay is approxim- 
ated as an all-pass filter [2] for which several parameters must be 
determined. This means that, for a simple system, the total number of 
polynomial coefficients which must be evaluated is greatly increased 
by the addition of a time delay. Because of this disadvantage, a 
different approach was adopted for fitting parameters to the motor 
unit transfer function estimates. 

DE HCE,) is the transfer function estimate and Gh), is the 
transfer function which is fitted to this estimate, a complex error 
term e, may be defined as 


k 


e, = AGE.) - HGE,) (6.8) 


oxsdw boxtmtata “ 
a! shar - 
ch | -S =* i}! 
eegysvn0s nokaestoxs elf? sont vodwun soksavett eit aneee: ay 
ak mist Youre ond dotw hegstooeus sotiomaygnisdetow ed +ye 03 wibiqet 
tanteo% bnu oenedjacs? .(00!) ywod yd beau eo (2t)9 anita gadte1 yitqu 
1O119 basdg bows ne 20 akdasd edt no baddtt aevans Jen? heey? enomeb 
atom médeve xolumo> 6 Yo stamiies notwonyt yelensrd off wollod mzr93 
at dotdw wed torte as 0) gaolbrosen bsgsli esvany ob mend yloaolo 
| C3109 datw besdgtew 
bra eios172 yd betiqqs nasd ead stubago0rq ayitave:! emf = e 
Hobie aredaye of e1osomeysq gakistt} te matecme ada od [OL] {Lirigmslt 


cl 


“mbxorqgn af etas soty sda .,dosoraqa afds daiw .velsb omky e sbulomt 
"iy 

ad daum ersionsisg Istevea doldw vet [S) 299059 aeseq-Ile as en base 4 
to ysdiwn isiod oft ,msiaye eslquie s yo! ,Jeda eamsom alal cbniieesniisd) 


bovesiont yfsnoxg ef baisulovs od teum dstdw aynenatdteca katwomedag ea 
B ,againevbasth elds to enysoat +¥etoh. amt 8 20 nots thm, Yd 
totam ot 04 arsiematagq ynlttit yo betgoba co doncsans aa 
.astania ee atreh a we 

oni et (ate bos. sdswtres notions valaneys edt a (apa i 
| Torre mabgnon © sea aida og vo at 


79 


The best-fitting transfer function over the frequency range f) ea 
m 


was defined as the transfer function for which € is minimized where 


(6.9) 


This quantity is a function of four parameters: k, 7, f , and t--the 

n 
time delay. Starting with a set of approximations for these parameters, 
e 


values for the fitted transfer function at frequencies f Piece 


ites geeks m 
were determined. The parameters were then modified according to the 
Fletcher-Powell method [110] until a minimum for ¢ was reached. The 
parameter values at this point were considered to be those of the 
transfer function which most closely approximated the transfer function 
estimate. 

The Fletcher-Powell method is a numerical technique for 
determining the minimum of a function of several variables. Because 
the minimum can be found quickly with this approach, the Fletcher- 
Powell method was selected in preference to other available minimiz- 
ation techniques [111]. Appendix 5 describes the details involved in 


applying this method to the problem of fitting transfer function 


estimates. 


6.2.5 Confidence Intervals for Fitted Parameters 

A computer simulation of the input-output relationship for a 
single motor unit was used to derive approximate probability distribut- 
ions for the fitted parameters. 

For the situation illustrated in Figure 6.4, the coherence 


between the input and output is [88] 
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2 oo Lvaeie a 
xy ~ 1+ K(f) (6.10) 


where K is the ratio of the power spectra of the uncorrelated noise 
signals, N,(£)/N, (£). The linear system was simulated by means of the 
digital second-order low-pass filter derived in Appendix 6 according to 
a method presented by Radar and Gold (112). With this approach, the 
filter parameters could be independently varied by adjusting the 
digital filter constants, and the coherence could be varied by changing 


the relative amplitude of the added noise, N The ratio N, (£)/N, (£) 


2° 
was constant at all frequencies, so the coherence of the simulated 
motor unit was independent of frequency. 

Approximate probability distributions for the various fitted 
parameters were obtained from 100 independent transfer function esti- 
mates for the simulated motor unit. These estimates were calculated 
according to the same procedure used to calculate the transfer function 
estimates for actual motor units. Specifically, records consisting of 
128 evenly-spaced samples of x(t) and y(t) were Fourier transformed, 
smoothed estimates were calculated by averaging the results from 
consecutive records, and estimates for the system transfer function 
and coherence were calculated from the smoothed spectral estimates. 

The parameters which best described these transfer function estimates 
were determined according to the technique described in Section 6.2.4, 
and histograms such as the one illustrated in Figure 6.5 were plotted 
to show the distributions of the parameters. The same procedure was 
repeated for transfer function estimates simulating a variety of motor 


units, and, in all cases, the shape of the histograms suggested that 


the fitted parameters are approximately normally distributed. This 


88) GineT * ex’ 
selon baxbLexzoom offi 3 jo /eitoege sewng. pil? vo obs add et > 01 
ei Yo ansem yd bogelomia cov maseye toot att. (i) NGI ge atargte 
64 gatbioogs a xbbasqqA ot bevitah reaflt seaq-wos sabte-bnossa fa3 bg tb 
of? .dosotdqs atta A2aW .ES00) bod bas tahal vd batnevetq bodiem a 
edi pnidvevtbs yd balrev yanabnegsbat ed bivoo stsiammxeq s93Llt 
Qnignads yd bolrev sd hitees sonsitoded 349 bap ,etostancs redL23 intigtb 
(2) pA) i otder odT agit 3 seten babba eda to shutiiqus evissisx od . 
bsiofomte aris Yo ssneTenos sit o8 “Bbroneleest tie Je Jneteqoo 8sw 
.yonsvpert Jo Jushnaqgsiat: esw tiou tosom 
beds? avotisy sf? sol andisudivieth yilitdedorg si amtxotgga 
-~tie9 notjoaqut tatensxt Snsbnaqsbat O0f most bontatdo eYaw s19JeqeTeq 
betaiuslas syaw astamitas susdT tine yoidn betsiustte ofa rod aetam 
coitoaut 19fems 13 9f1 sjalustso of baeu stubasoxg suse S23 OF gutivocan’ 
to gnisdeianos abto9s7 evileoittiosg? «vita x yoiom leytos t0F esiamttes 
~hewrolaasy) satrot araw (3)v bos (ad) To esfgose bsoagqe-_inevs est 
mort asloues1 adi pnksstavn vd bsisivolso s19W aajamiiee bedjooms , 
votiom? tslensys meteye sda 701 setemijes bas ,ebadse: svivuasencs, 
-eoiamistes Llaxiseqe bedtoome sft moti betsiuaiss syew sonezados bas 
sojaniies notjonnd tetene12 seadd badizoesb sau fotdy exsisuetsq odT 
»§.$.9 nokiss® at baditsasb supindost sf3 03 gatbropde bateeasian s19W 
batiolg staw <,.d srugitt Ar betstyeulit sno 303 as ded ansygotetd bas 
asw sYubs9otq ome odT .e1sJomeysg od? Io snotspdbrietb ons —_— os! 7 
toton to Kanssby & gotielumta estanites notiony? 1stenss3 oe. >s é 
dnd boseauaua ampagotaid si3 20 sqedé od3 ,eeen> fs at bas exo 
eh? .boandhazatb yLlamon Uldum ixoiqan 98m atvsamenag basa ai) 


4 


78 


4) 

see] 

em T 

RR KKK KH EEL 

2 Se He he OK a a ae HK He eH HH HO 
SKRSKKEKRKKKKK EKA HH KAKK KKK 
RRARKSSHKREEKSAHKKEKEARKH HK HK KREKEKSHE | 
SRKERHEEKAA REA EKAA ARSE RK RSE RAS ERK RAE EAH EHH KE HD | 
REKKKRKR EKA ER KEK KSA HK KHKK KH K K HK ] Y 
SRSEKSKHKAKKKREEKR SRT SKK RARE KKKARHEAE KARE K KE HT | 
KEKE KS SAH SKKHEKEKKKKKAK KS KRKKKKEKE KKK | 

ee SKK S HH HK KH HH HE HO HK HE 

kekegeekee eee ee RG 

eke Oe aL 

mk & | 

RRERKBS 


Soo 8 


ee eet a ee ea eo ee 
ats ; (om ] 


S8IUGIINIIQ JO Jaquinn 


0.8 


DC Gain 


FIGURE 6.5 Distribution of DC Gain 


: 
i ae 
~ wee > © a 


. sy 


79 


assumption was verified with a chi-squared goodness of fit test [14] 
for several randomly-selected histograms. The test indicated that the 
distributions described by the histograms could not be differentiated 
from normal distributions at confidence levels of 90% or higher. 

Figure 6.6 shows how the number of periodograms which are 
averaged to yield a single spectral estimate affects the distribution 
of fitted damping ratio for a system with a response similar to that 
of a fast motor unit. (The damping ratio for this system was 
-1.0, and the coherence was 0.8.) Results similar to this were obtained 
for all parameters and for systems with coherences ranging between 
0.5-1.0. It is apparent that a substantial decrease in standard 
deviation can be achieved by increasing the number of periodograms 
associated with each spectral estimate from 5 to 10, but the standard 
deviation decreases much less sharply when the number of periodograms 
is increased further. Because an increase in the number of periodograms 
requires a longer period of motor unit stimulation, spectral estimates 
based on the average of 10 periodograms represent a good compromise 
between the desire for low standard deviation and the need to keep the 
stimulation period as short as possible. Such estimates were used for 
all subsequent studies. 

Since the variances of spectral estimates for systems with 
low coherence are high, the parameters fitted to such spectral estim- 
ates have high standard deviations. This effect is illustrated for a 
typical situation in Figure 6.6 which shows that the standard deviation 
increases rapidly as the coherence decreases from 1.0. Expressions 
which show the effect of coherence on the distributions of each fitted 


parameter will be presented later in this chapter. 
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FIGURE 6.6 Effects of Record Length and Coherence on 
The Value of Fitted Damping Ratio 
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For single motor units, the squared difference between the 
fitted transfer function and the transfer function estimate was much 
smaller for high frequencies than for low frequencies. Since curves 
were fitted to the transfer function estimate so that the sum of these 
squared differences was minimized, the low frequency components of the 
transfer function estimate had a much greater effect on the sum than 
did the high frequency components. Therefore, when the high frequency 
portion of the meena ee function estimate was neglected during the 
fitting procedure, the effect on the fitted parameters was small. 
Because motor unit transfer function estimates tended to become unstable 
near the Nyquist frequency (as in Figure 6.2A) the last 20 points of the 
64 point transfer function estimate were not used in the fitting 
procedure. This maneuver did not greatly affect the distributions of 
the fitted parameters, and it yielded a reduction in the computing time 
required to complete the calculations. 

The scaling of statistical data so that it fits a standard 
distribution is a common procedure. An equivalent approach involves 
scaling the standard distribution so that it fits the data. Such a 
scaling procedure was applied to the normal distributions of fitted 
de gain and natural frequency. That is, standard distributions of 
fitted parameters were calculated for a simulated motor unit with a 
de gain of 1 and anatural frequency of 1 Hz. Distributions for motor 
units with any other values of k or im were obtained by scaling the 
standard distributions. This procedure was adopted when it was found 
that distributions obtained by scaling closely matched distributions 
which were obtained directly for a variety of simulated motor units. 


Distributions for fitted time delay were determined for 
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simulated motor units with delays equal to zero and half the sampling 
interval. Since the standard deviations of the fitted delay differed 
by only about 20% for the two situations, and since the time delays 

for actual motor units were usually about half the sampling interval, 
the distribution calculated for this value of delay was assumed to be 
valid for all responses of actual motor units. The distributions of 
the other fitted parameters were negligibly affected when the delay was 
changed from zero to half the sampling interval. 

Because the distributions for fitted damping ratio obtained 
by scaling did not match distributions which were obtained directly, 
the direct procedure was used to determine distributions for a number 
of different damping ratios within the range observed for actual motor 
Uni esc’ A nea tintereot geton technique, to be described below, was 
then used to calculate distributions for values of damping ratio other 
than these. 

The distributions of the fitted parameters depended upon the 
frequency range covered by the fitted transfer function. This effect 
is illustrated in Figures 6.7 and 6.8 for a simulated motor unit with 
the following parameter values: k=l, t=1, fl Hz, and T=25 msec. The 
Nyquist Pasecnay for the transfer function estimates was 10 Hz, and the 
upper limit of the frequency range spanned by the fitted transfer 
function is designated as Eve The plots show that moderate changes in 
the ratio f/f domot greatly aktecrt the distributions of the fitted 
parameters. 

For actual motor units, the ratio a ranged between about 
0.06-0.14. Distributions were determined for a simulated motor unit 


with a ratio of 0.06875, and these distributions were assumed to be 
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valid for all parameters fitted to the transfer function estimates of 
actual motor units. 

Also illustrated in Figures 6.7 and 6.8 are the effects of 
changing the frequency resolution of the transfer function estimate. 
Although all of the oe ae for actual motor units were determined 
with a frequency resolution of £,,/40 Hz, the figures show that, if 
necessary, reasonable results may be obtained with a frequency 
resolution as low as f,,/ 20 Hz. 

The plots in Figures 6.7 and 6.8 indicate that the relation- 
ships between bias and coherence and the relationships between coeff- 
icient of variation and coherence are approximately linear for 
coherences between 0.6 and 1.0. The coefficients describing these 
relationships were calculated by applying linear regression analysis 
[113] to data resulting from fits to different simulated motor units 
with coherences between 0.6 and 1.0. If n is the mean of a fitted 
parameter (k, f or ©), and € is its standard deviation, the following 
relationships apply: 

Reo= Vala Aa (6.11) 
2 =—be =p obey (Gs 2) 

When distributions for systems with damping ratios between 
0.8-1.5 were determined, it was found that the mean and standard 
deviations for fitted damping ratio were well described by the follow- 


ing relationships: 


tl 


nT) a. + ay + ab (6:43) 


0 
etc = by + by + bio (6.14) 


Values for the coefficients in these expressions are listed in 


Table 6.1 together with correlation coefficients. 
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When curves are fitted to experimental data, the true 
coherence and damping ratio are not known. Estimates of these quan- 
tities must therefore be used in Equations (6.11-6.14). Coherence 
estimates were calculated for the simulated systems described above, 
and the results of these calculations indicated that the average value 
of the coherence estimate over the lowest 30% of the frequency range 
spanned by the estimate is a reasonable measure of the true system 
coherence. The bias in this coherence measure is less than 3%, and the 
standard deviation is less than 0.03 for true coherence values between 
0.6 and 1.0. The only available estimate of the true damping ratio is 
the fitted damping ratio, and the results for the simulated motor units 
indicate that the bias in this quantity can be as high as 6% and the 
coefficient of variation can reach 20%. Means and standard deviations 
predicted by Equations (6.13) and (6.14) must therefore be viewed only 
as approximations. 

Equations (6.11) and (6.12) must be scaled if they are to be 
used for finding the distributions of k, f. and t for systems in which 
these parameters have values other than those indicated in Table 6.1. 
The coefficients for time delay are valid for transfer function esti- 
mates based on data collected at 50 msec intervals. If a different 
sampling interval is used, these coefficients must be multiplied by 
. A/(50 msec) where A is the sampling interval. The coefficients for 
de gain and natural frequency are valid for systems in which these par- 
ameters have values of unity. For systems other than these, the 
coefficients must be multiplied by the true parameter value. Of course, 
in experimental situations, the true parameters are not known, so 


scaling must be implemented with the fitted parameters. The use of 


| ‘eure sits ,preb Lad nim rsgxs oj bsioli sts qavaw5' ‘ane eatin 
-naup sasdd to eeipalied seni gon ors) bIRT ghtipend ns soastedos 
aopetsiiod =. @1.d-i1-o) “epolgnugh nm bagy sd srotereds J erm aotala 
,svods badiiysesh amsteye 4S sinmbe add. 164 nei a tse oraw essnmiies | 
auiay auntevn oft Jers batenion ansitelesian. sagt eo ativess oft bas 
ages yorelpst? ads 20° UE syawol silt 1avo, sJantis= sonexados ‘edd to 
isseVe Suto oft, to aziebon afdsnoesoy & at sismizes oda Ye bonne qe 

ads bor. $C nad ‘gael ‘ab styeesm sanszsies .ify ak anid s#% .sansitsioo 
qabwisd esuiey somstSo> Su33_ rot £0,0 cet eesf ef nottatveb byabasse 
ef okja: gniqesl sud wii, to siamgoes J folretes gio oft 40,1 bam 0.0 
Sitqu totum badaiunse adi iso etivgs of? boe ,obse1 anigmsb be3252 oft 

\ 

git bne SO en rigid ee’ oc n65 vitbasup aily ol entd.sft Jedd sisotbnt 
envigatves Disboese hrs enseM AOS “doser (59 sotisiasy to trekgttisos 
foo bswaty od stotsi sil) deun (Afea) bie (EL-8) enoigeups vd bestotbsag 


.snoitsemrnotqqs 8&6 


ad iod sts yous 22 GE - puum €S1.a) bew (12,8) amolioupa 
idw nl snedeve tot vol i to daotiudivsarb 9a, sivbbaks +07 beau 
sidel arid thant seody Gend sano asulev sveil erat aemeg sears 


aes noijnnud weienes) toi | ey aie vadeh- said Yor giasrstitess sft 
ineszotit! i -igpoejab sgsen DO? 9 bornefios Aish 16 bseed e@otea 

“at bal borslow gd Jem vinatsiti sos weary (boa et tavassat ime biaii 

101 f4neiciiteses 4d .levisiot atttignse gid eD 5 state deal 02)\4 . 

| -1sy seany ibadv nt emsteve 20% Srlsy sys eli ve Laake bag afeg ob 
of4 ,. oes? nad’ wolto. emadaye 10o% ,yriny, to Belay aved atstsms 


ssetioo 20 ~ .oulsev got MEISE * enix, Brit yu! fd teal ea ad Jeu aii 


oz ,nwond ton 9tB eros ezsg sur ‘al aed fosaamtvagee a 


¥o. se) SAT .erstemeyseqg beti1ii eda d3tw dedaemalene sd eu i Dee 


88 


such fitted values in scaling is the major source of error affecting 


the estimation of distributions for fitted parameters. This error is 


introduced because, for a system with a coherence of 0.6, the 


coefficientof wardation for fitted de gaimris: aboutrad 4%. tiFor fitted 


natural frequency, the coefficient of variation is about 9%. The bias 
in fitted dc gain and natural frequency is about 1% and 3% respectively. 
The coefficients in Table 6.1, when properly scaled, were 


used to calculate estimates of the means and standard deviations of 


the distributions for fitted parameters. Since these fitted parameters 


were assumed to be normally distributed, approximate confidence 


intervals for each fitted value could be calculated. Even though these 


confidence intervals were only approximate, they were useful in 


interpreting the transfer function estimates for single motor units. 
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CHAPTER 7 
RESULTS 

7.1 Classification of Motor Units 

Smith and Ovalle [34] classified muscle fibres in the 
iliofibularis muscle of Xenopus laevis into five groups according to 
histochemical and ultrastructural characteristics. Smith and Ovalle [34] 
also observed that some muscles do not contain all five fibre types 
found in the iliofibularis. Studies by other investigators [35,71] had 
indicated the existence of histochemically-distinct muscle fibre types 
similar to those identified by Smith and Ovalle [34]. As well, 
Lannergren and Smith [35] had demonstrated a relationship between the 
histochemical and functional properties of muscle fibres; and Smith and 
Lannergren [72] had classified motor units into four groups, two fast 
and two slow, according to functional characteristics. Smith and 
Lannergren [72] had also concluded that all the muscle fibres comprising 
a given fast motor unit are histochemically and functionally similar. 

Based on the above evidence, Smith and Ovalle [34] proposed 
that skeletal muscle fibres in the Anura may be separated into five 
structurally and functionally distinct classes consisting of three 
types of fast muscle fibre and two types of slow muscle fibre. An 
advantage of this hypothesis is that it can be tested experimentally. 
If it proves to be correct, the problem of understanding the organiza- 
tion of the motor control system will be greatly simplified. 

To test the hypothesis of Smith and Ovalle [34], an attempt 
was made to classify the motor units studied in this investigation into 
five groups, three fast and two slow. It was felt that if an approp- 


riate classification scheme could be developed, if it could be shown 
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that the groups are indeed different from one another, and if a corres- 
pondence between the motor unit types and muscle fibre types could be 
demonstrated, the hypothesis would be greatly strengthened. 

Several criteria were considered during the classification 
procedure. .- Most important, were qualitative differences among the 
responses of the motor units to various stimulus signals. Since 
specific muscle fibre types appear in specific areas within the muscle 
[34], the location of a motor unit within the muscle was an important 
factor. Also of interest were quantitative differences in particular 
characteristics of a motor unit's response. Finally, since Smith and 
Lannergren [72] had previously classified motor units into groups, 
their classification criteria were considered with some interest. 

The first step in the classification procedure involved the 
separation of motor units into fast and slow categories. Fast motor 
units were considered to be those which produced a twitch response when 
subjected to a single stimulus pulse. Those which responded only to 
repetitive stimulation were classed as slow motor units. 

All slow motor units produced a gradual increase in tension in 
response to regularly repetitive stimulation, and a gradual decrease in 
tension when stimulation ceased; but some motor units responded much 
more quickly than others. As well, if stimulation was maintained, 
tension produced by the units which responded more quickly reached a 
plateau and then began to rise again. If stimulation was stopped 
during this second increase in tension, and a series of single stimuli 
were then applied, small discrete contractions could be observed. 

No such behavior was apparent in the slowly-responding motor units. 


Based on these differences, slow motor units were classified into two 
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groups. This classification procedure for slow motor units was also 
used by Smith and Lannergren [72]. Because the type 5 muscle fibres 

of Smith and Ovalle [34] are believed responsible for the type of 
contraction exhibited by the more slowly contracting motor units [69], 
these motor ee were designated as type 5. As suggested by Smith and 
Ovalle [34], the remaining slow motor units were designated as type 4 
to match the type 4 muscle fibres. The initial tension ratio, defined 
in Chapter 6, was helpful in separating the slow motor units. Initial 
tension ratios for the type 5 units were all less than 0.25, whereas 
those for type 4 units were greater than 0.45. 

A type of response not observed by Smith and Lannergren [72] 
was apparent for several of the fast motor units studied. The twitch 
responses for these units were much smaller than for other fast motor 
units, and consecutive twitch responses varied widely in amplitude. 
These motor units were designated as type 3 because they were always 
located near the core of the muscle, and the only fast muscle fibres in 
this area were those designated as type 3 by Smith and Ovalle [34]. 
Although the twitches for type 3 motor units were small, the maximum 
tetanic tensions (defined in Chapter 6) often approached values similar 
to those of the other fast motor units. Therefore the twitch/tetanus 
ratios (defined in Chapter 6) for type 3 units were much smaller than 
they were for the other fast motor units. The twitch/tetanus ratios 
for type 3 motor units were always less than 0.09. 

The greatest differences among the remaining fast motor units 
were in the rates at which the units fatigued when subjected to repet- 
itive stimulation. Following Smith and Lannergren [72], these motor 


units were separated into two categories on the basis of their resis- 
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tance to fatigue. Since all these motor units were located away from 
the core of the muscle, and since the only muscle fibres in this area 
were the types 1 and 2 of Smith and Ovalle [34], all fast motor wits 
not classified as type 3 were classified as types 1 and 2. The more 
rapidly-fatiguing motor units tended to be those which produced the 
largest twitch responses, and these motor units were located in the 
outer portion of the muscle. They were designated as type 1 because 
type 1 muscle fibres were the largest, and on the basis of their histo- 
chemistry, the most susceptible to fatigue. Also, type 1 fibres were 
found only in the outer portion of the muscle [34]. The remaining fast 
motor units were designated as type 2. They were located nearer the 
centre of the muscle, as were the majority of type 2 muscle fibres. 

The fatigue coefficient (defined in Chapter 6) for type 1 motor units 
was less than 0.75. For type 2 units, this quantity was greater than 
O8S2. 

It should be stressed that the above classifications were 
only tentative. However, because the results of frequency response 
studies and histological investigations, to be reported later in this 
chapter, supported the classification scheme, no attempt was made to 
revise it. 

The various descriptive parameters defined in Chapter 6 
served to quantitatively describe the characteristics of the various 
motor unit types. The values of these parameters for the five motor 
my types are listed in Table 7.1 together with other characteristics 
which were observed. In addition, typical responses of the various 
motor unit types are illustrated in Figure 7.1. 


Type 1 motor units were supplied by the largest motor axons. 
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MAXIMUM TWITCH 
fare AXON TETANIC TENSION | AMPLITUDE | TWITCH- TWITCH FATIGUE 
tyrt | DIAMETER (vm) (g-force) (g-force) |TETANUS RATIO | PEAK TIME(ms) | COEFFICIENT? | \isgce oF 
TYPE | MEAN | S.D. MEAN [S.D. |MEAN [S.D. | MEAN OBSERVATIONS 
2 15.6 4.9 19.0] 11.7 2.7) 1.24 0.082 fo 12 
~ MAXIMUM INITIAL 
AXON TETANIC TENSION | TENSION 


MOTOR 


unit | DIAMETER (um) 


(mg force) 


RATIO NOTES: 


NUMBER OF 


TYPE [ MEAN 


OBSERVATIONS 


1. Based on 6, 6 and 4 


observations for motor 
unit types 1, 2 and 3 


respectively. 


TABLE, 7.1 


Characteristics of Single Motor Unit Twitch and Tetanus Responses 
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FIGURE 7.1 Responses of Single Motor Units to Stimulation 
with Trains of Evenly-Spaced Pulses 
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They developed high tetanic tensions, but they fatigued rapidly. Of 
all the fast motor units, the twitch/tetanus ratios for type 1 units 
were highest. 

Type 2 motor units developed smaller tetanic tensions than 
type 1 units, and they exhibited lower twitch/tetanus ratios. They 
were supplied by medium-size motor axons. With maintained stimulation, 
the decrease in tetanic tension due to fatigue occurred more slowly 
than for the type 1 units. 

Type 3 units were innervated by axons with diameters ranging 
from 10-15 wm. Their resistance to fatigue often approached that of 
type 2 motor units. The coefficients of variation for twitch amplitude 
were calculated for a series of fifty consecutive twitches of a type 3 
motor unit at various muscle lengths. The results of these calculations 
are plotted re. with the mean twitch amplitudes in Figure 7.2. 
This plot illustrates that the twitch amplitude became less variable as 
muscle length increased, and that the mean twitch amplitude was highest 
when the muscle was near its natural resting length, 1) 

In one experiment involving a type 3 motor unit, a small 
biphasic longitudinal current record was recorded at the surface of the 
muscle in association with the twitch response. On tetanic stimulation, 
the amplitude of the longitudinal current record increased as the 
tetanic tension developed. If this stimulation was continued, the 
muscle tension eventually began to drop, presumably because of fatigue, 
and the drop in muscle tension was accompanied by a decrease in the 
magnitude of the longitudinal current. This suggests that the fatigue 
effects may be largely electrical in origin, and that they are not the 


result of failure in the muscle contractile mechanism. A similar set 
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of measurements for a type 2 motor unit agreed with observations made 
by Lannergren and Smith [35]. These observations indicate that the 


contractile mechanism begins to fail before the electrical activity is 


noticeably affected. 


7.2 Frequency Response 


Transfer function estimates were calculated for several 
examples of each motor unit type. Typical examples of these transfer 
function estimates are plotted in Figure 7.3 for the slow motor unit 
types and in Figure 7.4 for the fast motor unit types. ee shown in 
Figure 7.4 are curves representing the best-fitting transfer functions 
calculated according to the method outlined in Chapter 6. The horizon- 
tal line on the coherence plots represents the average coherence over 
the frequency range spanned by the line. For the fast motor units, the 
coherence remained high for frequencies below about 15 Hz, and a second- 
order model seemed appropriate for this frequency range. Above 15 Hz, 
the drop in coherence could be attributed to the effects of quantization 
error as discussed in Chapter 6. Because the last 20 points of the 
transfer function estimates were neglected during the fitting procedure, 
this drop in coherence had little effect on the shapes of the fitted 
curves. Curves are not fitted in Figure 7.3 because the low coherence 
associated with the slow motor unit types indicated that the transfer 
functions for these motor units would not be appropriately represented 
by a linear model. 

The results of fits to transfer function estimates for the 
three fast motor unit types are listed in Table 7.2. These transfer 


function estimates were determined with the muscles at their natural 
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FIGURE 7.3 Transfer Function and Coherence Estimates 
for Slow Motor Units 
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FIGURE 7.4 ‘Transfer Function and Coherence Estimates 
for Fast Motor Units 


resting lengths. The interpulse intervals of the stimulus pulse train 
were approximately exponentially distributed with a mean of 0.1 sec. 
Interpulse interval distributions describable by gamma distributions 

of orders 1-8 were tested on a few preparations. These variations in 
the order of the gamma distribution seemed to have little effect on the 
transfer function estimates, so the exponential distribution was used 
for all subsequent studies. Stimulation and sampling were controlled 
so that bursts of stimulation 2.56 sec in length occurred every 20 sec. 


Comparison of transfer function estimates calculated from the 
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first and last halves of groups of 20 power spectrum estimates obtained 
with continuous stimulation and with stimulus bursts proved to be an 
inconclusive means of demonstrating the advantages of stimulus bursts 
in reducing the effects of fatigue. This comparison indicated that the 
transfer function estimates were insensitive to such differences in the 
stimulation procedure over the time required to collect 20 records. 
Therefore, the justification for using stimulus bursts is based solely 
on the subjective observation that, with this technique, more samples 
could be collected over the course of the experiment before the reduc- 
tion in twitch amplitude due to fatigue became unacceptable. Generally, 
a 20% decrease in twitch amplitude over the course of the experiment 
was considered acceptable 

The values for each parameter fitted to the transfer function 
estimates were grouped according to motor unit type, and analysis of 
variance [138] was used to test the null hypothesis that, for each 
parameter, the groups came from the same population. This procedure 
uses an F test [138] to compare two quantities. The first quantity 
gives the average dispersion of the items in each group about the group 
means and the second quantity gives the dispersion of the group means 
about the mean of the entire population. If the null hypothesis is 
correct, these two quantities should not be significantly different. 
For all four fitted parameters, the F test indicated a difference at 
the 99% level of significance, thereby indicating that, for each 
parameter, at least some of the group means were different from one 
another. To discover which means were different from one another, 
comparisons between pairs of fast motor unit types were made with 


the least significant range test [138]. At the 95% level of signifi- 
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cance, the means of all parameters for type 3 motor units differed 
from those of the corresponding type 1 motor unit parameters. 
Differences at this level of significance were also apparent for all 
parameters except dc gain and time delay when means for type 2 and 
type 3 motor units were compared. A comparison of means for type l 
and type 2 motor units indicated that the natural frequencies and time 
delays differed at the 95% level of significance. 

The effects of changes in the muscle length and stimulus 
rate are illustrated in Figures 7.5-7.8 for examples of each fast 
motor unit type. Each figure consists of four separate plots. The 
two upper plots show the fitted parameter values joined with straight 
lines for each motor unit type, and the two lower plots indicate the 
approximate 95% confidence regions for the fitted values. These 
confidence regions are the areas between sets of straight lines 
joining respectively the upper and lower confidence limits for each 
set of data points. The confidence limits were calculated according to 
the procedure described in Chapter 6. 

The effects of changes in muscle length were studied for a 
stimulus signal having a mean stimulus rate of 10/sec. Transfer 
function estimates were obtained for a range of muscle lengths which 
roughly corresponded to the range of possible lengths for an iliofibu- 
laris muscle in vtvo. 

With the muscle at its natural resting length, transfer 
function estimates were obtained for a variety of stimulus signals with 
mean stimulus rates ranging between 0.2-20/sec. 

As Figure 7.5 indicates, the confidence regions for time de= 


lay were so wide that a systematic variation of this parameter in 
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response to changes in length or stimulus rate could not be detected. 
Values for fitted time delay ranged all the way from 3 to 18 msec, but 
in most cases, the fitted delay was in the 10 msec range. For the 
three preparations in which delay was measured directly, values of 

8, 10, and 11 msec were obtained. 

The delay consists of three components: the time required to 
transmit an action potential along the nerve from the point of stimul- 
ation to the neuromuscular junction, the transmission delay at the 
neuromuscular junction, and the delay associated with excitation and 
contraction in the muscle. For axon diameters between 5 Um and 20 um, 
the value of the first component ranged between 1-4 msec [116], and 
the synaptic delay was 0.5-0.8 msec [117], so the value of the 
remaining component had to be near 10 msec. An experimentally- 
determined value in the 10 msec range for this last component of the 
time delay is apparent in the results of Rudel and Taylor [46]. 
Differences among the mean time delays for the various motor unit types 
were too great to be accounted for completely by differences in the 
travelling times of the nerve action potentials. It therefore seems 
likely that differences in excitation-contraction delay exist among 
the various muscle fibre types. 

Figure 7.6 shows that the de gain for the type 1 motor unit 
was highest at very low stimulus rates whereas the de gain for the type 
2 and type 3 units reached a peak near 5 stimuli/sec. The de gain for 
the type 3 motor unit was highest when the muscle was near its nace al 
resting length, whereas the gain for the type 1 and type 2 motor units 
continued to increase as the muscle was lengthened. 


The damping ratios for all three fast motor unit types, 
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plotted in Figure 7.7, increased slightly as the stimulation rate was 
increased, and the damping ratios for the type 1 and type 2 motor units 
increased slightly for increasing muscle length. For the type 3 motor 
unit, the damping ratio was very high when the muscle was short, and it 
decreased sharply as muscle length was increased. When the muscle was 
2 mm longer than its resting length, the damping ratios for all three 
motor unit types were near 1.25. 

Figure 7.8 shows that the natural frequencies for all three 
fast motor unit types were highest for short muscle lengths or low 
stimulus rates. The value of this parameter decreased as muscle length 
or stimulus rate was increased. 

Values of the coherence estimates for each motor unit type 
are plotted in Figure 7.9 together with approximate 95% confidence 
intervals derived from the computer simulation of the fitting procedure. 
The plotted values represent the average of the coherence estimate over 
a frequency range between about 0.4 Hz and 12 Hz. The coherence 
for the type 1 motor unit remained above 0.8 for all muscle lengths 
studied, and the coherence for the type 2 motor unit remained near this 
value for muscle lengths below 1 ot mm. At the longest muscle length, 
the coherence for the type 2 motor unit dropped to near 0.6. For the 
type 3 unit, the coherence was less than 0.5 when the muscle length 
was short. As muscle length was increased, the coherence value 
gradually increased to a plateau of about 0.8 for a muscle length of 
1 +1 mm. The coherence for all three motor unit types in response to 
widely-spaced single stimuli was near 1.0, but this value dropped 


steadily as stimulus rate was increased. 
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Note: The vertical bars indicate 95% confidence 
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FIGURE 7.9 Effects of Mean Stimulus Rate and 
Muscle Length on Average Coherence 


. -¢ 


(>92\ mite) oto 2ulumit2 nosh 


ee ceccieimenncemarty t itr ~ RO 


snnsbi Laos 8ee sanatbat ard [ooroasy, ait 
1&0 bre. VO) Ses Be ‘bo ASS TSHIS yaoi alevaejni 


oe = a 


7.3 Mechanical Properties 


The transfer functions relating small changes in muscle 
length to changes in tension were similar to those obtained by Halpern 
and Alpert [104] for the sartorius muscle. A typical example, illus- 
trating the response to random length changes never greater than about 
80 um, is shown in Figure 7.10. This maximum displacement is similar 
to that observed for the so-called "isometric" force transducer during 
the measurement of motor unit transfer function estimates. For freq- 
uencies below 20 Hz, the coherence was near 1.0, the gain was constant 
to within about 2 db, and the phase changed by less than 20°, This 
means that the mechanical response of the muscle can be approximated 
by that of a spring for frequencies below 20 Hz. Each transfer func- 
tion estimate which was calculated for a single motor unit can there- 
fore be viewed as the response of the motor unit's muscle fibres in 
parallel with a spring representing the remainder of the aan At 
low frequencies, the compliance of the whole muscle ranged between 
about 0.1-0.35 mm/g-force for muscle lengths between 1-2 mm and 
1 +2 mm. 

| iG 1/K, is the compliance of the force transducer, 1/K, is 
the compliance of the muscle, and F(t) is the force generated by the 
single motor unit, the force indicated by the transducer will be 
Pie ute Git 
K TK, 
Therefore, for a force transducer compliance of 0.01 mm/g-force, the 
maximum attenuation of the force generated by any motor unit was only 
about 7% and this attenuation was constant for all frequencies below 


20 Hz. For this reason the effect of the mechanical properties of the 
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FIGURE 7.10 Relationship Between Muscle Length and Tension 
for Small Random Perturbations in Length 
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passive muscle fibres on the frequency response of single motor units 


was not subtracted from the transfer function estimates. 


7.4 Histochemical Studies 

The photographs in Figure 7.11 show the results of studies 
in which single fast motor units were continually stimulated in order 
to deplete the glycogen in their muscle fibres. These studies were 
carried out on 2, 3, and 3 examples of motor units of types 1, 2, and 
3 respectively. Certain muscle fibres in the photographs of sections 


Stained for SDHase are marked with numbers which indicate the fibre 


type according to the designation scheme used by Smith and Ovalle [34]. 


The same fibres in the sections stained for glycogen are marked with 
letters indicating the degree of glycogen depletion. The letter U 
indicates a fibre in which no depletion is apparent, P indicates a 
partially depleted fibre, and D indicates a fibre for which the degree 
of depletion is substantial. Depleted, partially depleted, and undep- 
leted fibres of the same type were found in each preparation, but the 
numerical designation for those fibres exhibiting any evidence of 
glycogen depletion always matched the numerical designation applied to 
the isolated motor unit. No depleted fibres were found in a control 


preparation which was not stimulated before being fixed. 
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FIGURE 


Vere 


Transverse Sections of Muscle Fibres From Three 
Different Nerve-Muscle Preparations. 


The isolated motor units were identified as 
Type 1 (A,B), Type 2 (C,D), and Type 3 (E,F). 
Serial sections were stained for SDHase 
activity GA,G,E) and glycogen (8;D,F).. The 
scale bars indicate 200 um, 
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CHAPTER 8 
DISCUSSION 
8.1 Significance of Transfer Function Estimates 

This investigation was designed to determine if spectral 
analysis methods can be profitably applied in studying the input-output 
relationships of single motor units. The iliofibularis muscle of 
Xenopus laevis was known to contain motor units of several different 
types, and it was felt that the effectiveness of spectral analysis 
could be assessed from the results of in vitro studies involving these 
motor units. It was also hoped that such studies would provide infor- 
mation concerning the function and organization of the various motor 
unit types. 

Spectral analysis was used for calculating transfer function 
estimates from which linear models could be derived. A simpler 
approach to estimating the transfer function would involve Fourier 
transforming the waveform of a single twitch. This approach is 
possible because a single twitch is equivalent to the impulse response 
of the system, and the impulse response is the inverse Fourier trans- 
form of the transfer function [118]. Unfortunately, as the results in 
Chapter 7 indicate, such a transfer function estimate is a poor 
representation of the response to stimulation at higher rates. There- 
fore, the use of spectral analysis to study the responses for a variety 
of stimulus rates seems reasonable. 

The linear models presented here for the stimulus-response 
relationships of single motor units have several limitations. They 
are only applicable to certain motor unit types, they are useful only 


for predicting isometric muscle tension, they depend upon such factors 


as: 
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as stimulus rate and muscle length, and they neglect the effects of 
fatigue. Despite these limitations, however, such models are useful 


because they provide a standard basis for describing experimentally- 


observed responses. The parameters associated with each model can be 


used to compare the response of a particular motor unit with the 
responses of other motor units, with theoretical responses derived from 
mathematical models, or with responses of the same motor unit under 
different experimental conditions. It should be stressed that the 
linear models presented here are valid only for signals having 
perturbations in interpulse interval similar to those of the stimulus 
signals used in developing the model. For such signals, the standard 
deviation of the interpulse interval is equal to the mean interpulse 


interval (135). 


8.2 The Coherence Estimate 

Coherence estimates for fast motor units differ from those of 
linear, noise-free systems in two ways: the coherence drops sharply 
at high frequencies, and, although the coherence at low frequencies is 
approximately constant, its value is less than unity. 

Quantization noise seems to be the major cause of the drop in 
coherence at high frequencies, and this means that the coherence 
estimate could be improved by pre-whitening the system output signal 
[73]. Pre-whitening involves passing the signal through a well- 
designed filter so that the spectrum of the resulting filtered signal 
is approximately flat over the band of frequencies below the Nyquist 


frequency. The level of this signal may then be adjusted so that 
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components at all frequencies are digitized using the maximum resolut- 
ion of the analog-to-digital converter. Of course, the effects of the 
pre-whitening filter would have to be removed from the resulting 
transfer function estimate. Pre-whitening is not required for the 
input signal because a signal with a flat spectrum may be chosen as an 
input. 

The results in Chapter 7 suggest that the coherence for a 
fast motor unit is approximately constant with a value close to that of 
the low frequency coherence estimate. This value is less than unity 
presumably because of the combined effects of added noise and system 
nonlinearity. The response of type 1 and type 2 motor units to trains 
of uniformly-spaced impulses did not appear to be contaminated with 
added noise. No spurious or missing twitches were apparent, and the 
twitches were all the same size. For higher stimulus rates resulting : 
in fused muscular contractions, no sudden changes in muscle tension 
could be observed, Similar observations applied for type 3 motor units 
at long muscle lengths. For short muscle lengths, however, the coeff- 
icient of variation for twitch amplitude was high, and tetanic muscular 
contractions were not smooth. The low frequency coherence for the type 
1 and type 2 motor units was not greatly affected by changes in muscle 
length, but, for the type 3 unit, this value dropped sharply at short 
muscle lengths. It therefore appears that the coherence of fast motor 
units. is less than unity because of the nonlinear response exhibited 
by these motor units. In addition, the coherence of type 3 motor units 
at short muscle lengths is further reduced because of the noise-like 
effects which are observed under these conditions. 


Experimental measurements of the responses of a few examples 
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of each motor unit type have led to the assumption that the response of 


a fast motor unit is similar to that of a second-order low-pass filter 
with a time delay. Based on this assumption, simulation studies of the 
type described in Chapter 6 have provided useful information concerning 
the application of spectral analysis to the study of fast motor units. 
Such simulation studies would not be needed if experimental data could 
be obtained for a large number of single motor unit preparations, but 
such preparations are difficult, and the amount of information which 
can be obtained from each one is limited because of the effects of 
fatigue. 

The responses of simulated motor units only crudely approxi- 
mate the responses of real motor units, but these crude approximations 
are sufficient for studying the spectral analysis procedure. 

Neglecting the effects of fatigue, the greatest difference between the 
responses of the real and simulated aes units can be related to the 
coherence estimate. For real motor units, the coherence is reduced 
from unity largely because of the nonlinear relationship between 
stimulus rate and muscle tension. For the simulated motor units, the 


coherence depends upon random noise which is purposely introduced into 


an otherwise linear stimulus-response relationship. However, the 
probability distributions of transfer function estimates for a noisy 
system and a nonlinear system are identical if the coherences of the 
two systems are the same [23] (so long as the same spectral analysis 
procedure is used in both cases). Therefore, parameters fitted to a 
transfer function estimate for a simulated motor unit will be statis- 
tically equivalent to parameters fitted to a transfer function estimate 


for a real motor unit with the same coherence. 
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8.3 Relationships Between Motor Units and Muscle Fibres 

Because the location of each motor unit type within the 
muscle matches the location of a particular type of muscle fibre, and 
because the isometric contractile properties of the various motor unit 
types are unique, it appears that a motor unit contains muscle fibres 
of only one type. The numerical designation which was assigned to each 
motor unit type indicates the type of muscle fibre from which the motor 
unit is presumably composed [34]. Direct evidence for such a relation- 
ship was provided by the glycogen depletion studies. These direct 
results are based on a limited number of experiments, and not all motor 
unit types were studied; but, in every case, the hypothesis that single 
motor units are composed only of muscle fibres of one particular type 
was supported. Both motor units and muscle fibres can be classified 
into two major groups--fast and slow, but as more criteria are consid- 
ered, each of these groups can be subdivided. The possibility of 
further subdivision still exists, but, in the absence of any evidence 
suggesting the existence of additional sub-groups, the present classifi- 
cation of motor units into five distinct groups must be accepted. 

Although amphibian and mammalian muscle fibres are classified 
differently, the results of Weurker et al. [119] are applicable to this 
discussion. They proposed that single motor units in the medial 
gastrocnemius muscle of the cat are composed of muscle fibres which 
resemble one another very closely, and they suggested that the 
similarity among muscle fibres of the same motor unit is due to some 
type of neural influence. The results of cross-innervation experiments 


by Buller et al. [120,121], for the cat, and.tlulger al.wld22415 for-the 


frog, support this theory. 
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Other similarities between mammalian and amphibian motor unit 


organization are apparent. The largest motor units are, in general, 


innervated by the thickest nerve fibres. This relationship was 
observed in the cat by McPhedran et al. [57] and by Weurker et al. 
[119]. Since, in the froe at least, the largest motor units tend to be 
those which contract most rapidly, and since the conduction velocity of 
a nerve increases with its diameter [116], the relationship between 
motor unit size and axon diameter insures that the most rapidly 
contracting muscle fibres are capable of receiving stimuli as quickly 
as possible. 

In both mammals [119] and amphibians, the largest, fastest 
motor units fatigue more rapidly than do the smaller motor units. The 
large differences in the fatigue rates of fast muscle fibres in the 
frog, noted by Grabowski et al. [96], can be explained by the existence 
of three fast muscle fibre types with different fatigue rates, 
Lannergren and Smith [35] showed that fatigue in type 1 muscle fibres 
is not due to failure of propagation of the muscle action potential, 
and a similar situation is apparent for type 2 motor units. Fatigue in. 
type 3 units, on the other hand, develops very slowly and is accompan- 
ied by a decrease in the magnitude of the propagated action potential. 
It therefore seems possible that fatigue effects related to the 
propagation of the muscle action potential occur at the same rate for 
all fast motor units. The observed differences in the rates of fatigue 
for the three motor unit types are probably related to differences in 
the structure of their muscle fibres. The rapidly fatiguing type 1 
muscle fibres are poorly supplied with mitochondria, and their content 


of glycogen and fat droplets is low, whereas type 2 muscle fibres, 
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which fatigue slowly, are well equipped with all of these elements. 
Fatigue in the contractile mechanism of type 3 motor units occurs very 


slowly, and type 3 muscle fibres carry an extremely abundant supply of 


these energy-producing elements. 


8.4 Interpretation of Transfer Function Estimates 

The low coherences associated with slow noe units indicate 
that the responses of such motor units to random stimulation are poorly 
represented by linear transfer functions. This result is not unexpected. 
Earlier investigations [50,58,72] have shown that slow motor units with 
responses similar to those of type 5 units exhibit nonlinear character- 
istics such as the absence of a twitch response and the ability to 
maintain relatively high tensions at low stimulus rates. As Samosudova 
et al. [123] suggest, this nonlinear behavior probably results from the ' 
absence of propagated action potentials and the existence of a sparse 
sarcoplasmic reticulum which neither supplies nor reabsorbs calcium ions 
quickly. It is interesting that type 4 motor units also exhibit non- 
linear characteristics even though these motor units can shorten much 
more rapidly than can type 5 units. Although nothing is known about 
the function of single type 4 muscle fibres, one can predict that these 
fibres will be found to share many of the physiological properties of 
the type 5 Weeele fibre. 

In studying the function of slow motor units, the response to 
stimulation of a single nerve fibre is of interest; but, because these 
motor units are innervated by several motoneurons, a more detailed study 
of their function will require a different stimulation procedure. 


When muscle fibres contract in response to stimulation of 
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their motoneurons, a chain of complex processes is involved. However, 
for fast motor units, the combined effects of all these processes may 
be approximated with a second-order system and a time delay. The time 
delay appears to be associated with the events leading up to the 
release of calcium in the muscle and its subsequent contraction; and 
the second-order behavior can be simulated with the Huxley-Julian model 
for muscle activation and contraction outlined in Appendix 7. Although 
the model is a complex one involving several interacting variables, 
Stein and Wong [124] have shown that its frequency response is similar 
to that of a critically-damped second-order system. They proposed 

that this behavior is due to the limiting effects of two parameters. 
According to this proposal, the relaxation of the twitch is limited 

by the slower rate constant of activation (8 in Appendix 7), and the 
rising phase of the twitch is limited by the rate constant for forming 
cross bridges (f,)5 the rate constant for breaking cross bridges (g,) 
or the series stiffness (k_). 

Since the input signals to motor units normally consist of 
trains of narrow pulses, a study of the impulse response of a second- 
order system can be helpful in interpreting the frequency response 
plots which were measured experimentally. The form of the impulse 
response depends upon whether the damping ratio is less than or 
greater than unity [2]. Plots of the impulse response, hit) 4) non 
various values of t are illustrated in Figure 8.1. The response for 
G=1 is of particular interest. For this situation, known as critical 
damping, the impulse response decays from its peak value as rapidly as 


possible while still remaining monotonic. 


The effects of k and f on the impulse response can be 
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determined from Figure 8.1. An increase in k will result only in an 
increase in the overall Magnitude of the impulse response, thereby 
increasing the average value of the response. An increase in f 
causes the peak value to increase and the time course of the response 
to shorten, but the average value remains unchanged. 

Certain characteristics of the response to a train of 
narrow pulses can be deduced from the shape of the impulse response. 
If the natural frequency is low, summation of the individual impulse 
ee igeisers will be apparent at low stimulus rates, and the output wave- 
form will be relatively smooth. For such a situation, the output at 
any time will be noticeably influenced by impulses which occurred 
earlier in time, and the response to a sudden change in stimulus rate 
will therefore be slow. If the natural frequency is high, however, the 
output waveform for low stimulus rates will appear more like a series 
of distinct impulse responses, and the response to a sudden change 
in stimulus rate will occur more quickly. 

For all ae the longest muscle lengths, the responses of 
type 1 and type 2 motor units were close to being critically damped. 
The slightly higher damping ratios which occurred at long muscle 
lengths and for type 3 motor units are indicative of a more slowly 
declining twitch which results in a slightly lower speed of response 
to changes in stimulus rate. 

The increase in damping ratio and the decrease in natural 
frequency which occurred for type 1 and type 2 motor units as the 
muscle was stretched indicate that, as the muscle becomes longer, its 
speed of response decreases; and the higher dc gains at long muscle 


lengths show that the average force generated per stimulus pulse is 
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higher. These observations agree with the results presented by Close 


[94] for the frog sartorius muscle which contains only type 1 and 


type 2 muscle fibres [34]. He observed that the maximum twitch tension 


and the twitch contraction time increased as the muscle was stretched, 
and he concluded that these increases are due to length-dependent 
changes in the contractile process. 

The unusual behavior of type 3 motor units at short muscle 
lengths (low coherence, high damping ratio, variability of twitch 
amplitude) may be due to stretch sensitivity of the motor nerve endings. 
The effects of such stretch sensitivity have been noted by Kuffler and 
Williams [62] and by Ralston and Libet [64]. Orkand [65] proposed that 
these effects are due to the small amount of neuromuscular transmitter 
substance released at the motor endplate by an individual nerve impulse. 
Hutter and Trautwein [125] have demonstrated that neuromuscular trans- 
mission is facilitated when the muscle is stretched. Such facilitation 
is apparent as an increase in the size of the local depolarization 
which occurs in the muscle fibre near the neuromuscular junction when a 
stimulus pulse is applied to the nerve. At short muscle lengths, the 
local depolarization may not be great enough to cause a propagated 
muscle action potential. However, as Nasledov [63] suggests, the local 
depolarization could produce a small local contraction of the muscle 
fibre, and the effects of such a local contraction may possibly be 
great enough to increase the damping ratio. The muscle fibres affected 
in this way would not contribute to the twitch response of the motor 
unit, so this response would be smaller than usual. Since the response 
of type 3 motor units at short muscle lengths is highly unpredictable, 


it seems probable that these motor units are normally activated only 


-_ . 


_ ae ae ; 

a “7 2° 

ryt ve Ole ae - 
eid sim te aN . 
com 


sosbrhyal=doghet ri Sub aye pene ole 
_ ci vy: 


sijann iiaie IW wd ru aojom & velit solvaded Eowewar ; | fat 
tisaiws to “alitdetiav ,.ol384 gaiqaad dgtd dy td ; dora aks, argue of 
eguibos sviec soto siz to whiv ideas legit 0? sub od Petes _ : 
bea veitaux va badom mead evan Pane te dovedds dome to ‘woeite adt _ 
isd bexogorg (2d] Geet .[d@) sadka bas 09 aes ed ban [88] adhere | 
+533 tuene ts sslustomerwen to Tovoms Llane, ald oF sub ste psoe3t9 98 suedi 
,satvem) svren Tuubtythab as yd etalahha tojon odd ae baeseles pate ery | 
-#n612 iwwiluoeumotwsr tht betcaienoneb sved [eSh] sYowiuay? ban 2ed3uH” | 
aolsetitiowy? dowa .bardoisvda af sloeum old nate bsond bitoe? et nokeekm 
cotisstraloqeb lace! sd Yo seta ant ot senstant an inezeqqa er. 


s nodw holjoot 2eludaumotvan vd tesn o1dlt ofoanm oft ot apace dotdw 


Ah lant SS | wal 
afta ,adigne!l atoeun ai0le JA . ayvven of?.o7 baliqgs al salug sutuntie mit 
iad Ae wey 
botsgeqesg 5 sevso ot dquohs tasty od Jom Yen avtisstisfogsb Iscol 
; ing 


frocl ait , 2teegaue [€a) vobsluait «o Tia bacon wee oo an 
4 AD sal - oa 


TosaM add, Io sudegeer dothws ori 08 ere 08 b 
a Mt 


senogoos sain Acar > fauan 7 aa! fone ed > 
einetiere 2 or. otoeam 4 
Bi a baa 


au Peat nd hog yor 


at muscle lengths greater than the natural resting length. 

Because the input-output relationships for fast motor units 
appeared to be noise-free when the muscle was at its natural resting 
length, the decrease in coherence observed for all three fast motor 
unit types as the mean stimulus rate was increased can be interpreted 
as a decrease in linearity resulting from the development of tetanic 
muscular contractions. Also common to all fast motor units at 
increased stimulus rates was a systematic decrease in natural frequency. 
When considered together with the increase in damping ratio which was 
only slight, this drop in natural frequency indicates Te the motor 
units respond more quickly to a single stimulus than to variations in 
some mean stimulus rate. 

The peak in the dc gain versus mean stimulus rate curve 
occurred at very low stimulus rates for the type 1 motor units and at 
higher mean rates for the type 2 and type 3 units. Since the dc gain 
can be viewed as a measure of the average force generated per stimulus 
pulse, it appears that type 1 motor units respond most effectively when 
they are stimulated at low rates, and type 2 and type 3 motor units are 
most effective when stimulated at higher rates. In other words, the 
type 2 and type 3 motor units seem to be best suited to the detection 
of variations in some mean stimulus rate, whereas the type 1 units 
operate most effectively in an on-or-off manner. One effect which 
might be associated with this difference in rate sensitivity is the 
facilitation which has been observed at the frog neuromuscular junction 
in response to repetitive stimulation [126,127]. Repetitive stimulation 
also affects the Rectal and contractile processes in frog muscle 


[128-130]. Stein and Wong [124] demonstrated with a mathematical model 
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for muscle net eenion and contraction that saturation of the calcium 
re-uptake process could be responsible for the observed changes in 
muscle frequency response at high stimulus rates. 

In summary, an overall view of the possible functions of the 
various motor unit types can be presented. Slow motor units seem to 
have little influence on the control of normal movement. They produce 
very little force, and they cannot respond quickly enough to influence 
movements which occur at normal speed. Probably, the major function of 
these motor units is the maintenance of posture. Type 1 units appear 
to be designed to produce large, rapid muscular contractions in response 
to short bursts of stimuli. These motor units cannot be used effect- 
ively for repetitive movements because they fatigue rapidly, so they 
are probably used only when a quick powerful movement is required. 
Repetitive movements such as those required for normal swimming are 
most likely produced by type 2 and type 3 motor units with the type 2 
units dominating at shorter muscle lengths. The mechanism for generat- 
ing an increase in muscle tension therefore seems to depend upon the 
type of movement involved. 


In human muscle, when force is gradually increased from zero, 


the smaller motor units are the first to become active [131]. At low 
levels of force, recruitment of additional motor units is the major 
mechanism for increasing force, but at higher levels, increased firing 
rate becomes more important [132]. For twitch contractions, however, a 
different mode of activation is apparent [133]. Twitch contractions 
involving only the large motor units are commonly observed. 

If frog motor units are activated in a similar manner, the 


tension developed during rapid, powerful movements must be controlled 
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mainly by the recruitment of additional type 1 motor units. Stereo- 
typed, repetitive movements may be slightly influenced by variations in 
the stimulus rate, but movements of this type probably depend largely 
on the recruitment of additional motor units when increased force is 
needed. Movements requiring fine control of force probably use both 


recruitment and variation of firing rate depending on the level of 


force being produced. 


8.5 Conclusions and Recommendations 

Several conclusions resulting from this investigation have 
been presented in Chapters 6, 7, and 8. These conclusions are 
summarized below: 

1. Motor units in the iliofibularis muscle of Xenopus laevts 
can be classified into 5 types--3 pote td 2 slow. Each , 
motor unit contains muscle fibres of one type only. The 
histochemical characteristics of each muscle fibre type are 
consistent with the observed behavior of the corresponding 
motor unit type. 

2. If a single motor unit is stimulated with a train of narrow, 
randomly-occurring pulses, spectral analysis methods can be 
used to determine a transfer function estimate for the motor 
unit. This estimate relates the low frequency components of 
the stimulus train to the resulting isometric muscle tension. 

3. The results of such spectral analysis studies indicate that 
slow motor units are poorly-modeled by linear, noise-free 
systems at frequencies associated with normal movement. 


4. Because the transfer function estimates for fast motor unLts 
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are similar to the transfer functions of second-order low- 


pass filters, a curve representing the response of such a 
filter with an added time delay can be fitted to each 
transfer function estimate. In this way, each input-output 
relationship for a fast motor unit can be represented by a 
linear model with four parameters. 

5. Because the input-output relationship is affected by such 
factors as muscle length and mean stimulus rate, a set of 
linear models is required to represent the motor unit's 
response to changes in these factors. Such linear models are 
only approximate, and they are subject to many restrictions. 
However, they are useful because they provide a standard 
basis for describing motor unit responses. They also provide 
information which can be used to supplement the results of 
studies involving stimulation with trains of regularly- 
occurring pulses. 

6. Computer simulations of fast motor units, based on the 
experimentally-determined linear models, can be used to 
demonstrate the effects of various modifications to the 
spectral analysis procedure. Such simulations can also be 
used to estimate confidence intervals for the parameters 


describing the linear models. 


Future applications of spectral analysis to the study of 
muscle function should yield useful quantitative information on the 
responses of fast motor units to various combinations of length, load, 


and stimulus pattern. As well, the effects of other factors such as 
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temperatu i i 
Pp ture or drugs could be quantitatively assessed with this technique. 
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APPENDIX 1 
LOW-PASS FILTER 


The circuit in Figure Al.1 consists of a voltage follower and 
two identical second-order low-pass filters cascaded together. The 


transfer function of the first filter is 


0 3 1 
H(s) = —— - (Al.1) 


R 
1 atic. pe, @-> = 2 25) c4e Bee a7 


i Dig R, 
3 R, Ry Ry aa ge aes 


If the component values are selected so that the filter frequency 


response is maximally flat, the transfer function has the following 


form 


H(s) = ——z+—___.,,- (Al. 2) 
1 _ + ae 
mE. (2n£ .) 


where a is the gain of the filter and f. is its cutoff frequency. The 
component values required for cutoff frequencies of 3, 15 and 30 Hz are 
fisted in Table All: 

When two such filters are cascaded, the resulting magnitude 


and phase characteristics are respectively 


M(f) = ' i 
Peg 
Jit Cade) 
P(f) =- cane | 
1- am) 
Cc 
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APPENDIX 2 
GENERATION OF RANDOM PULSE TRAINS 

Random pulse trains with interpulse intervals distributed as 
gamma distributions were generated by a Hewlett-Packard 2100 computer 
with a digital-to-analog converter. The interpulse intervals for each 
pulse train were derived from uniformly distributed pseudo-random num- 
bers created BF the power residue méthod [135]. The particular set of 
pseudo-random numbers created consisted of integers with values between 
1 and 32767. A series of numbers distributed as a gamma distribution 
could be derived from these integers in the following manner: 


Te Xs is any number (not necessarily integer) from a 


population distributed uniformly between x and “x” })"and 
max min 
1 *i*nin ( . ) 
ee we pe he Pee A2,.1 
V4 68 wS -xX_, ) 
max “min 
the probability density function for y is [135] 
= Be PY (A2.2) 
pW) = Be 
Now if 
= (A223) 
a yy + Yo ate ane Tg 
the probability density function for t is [135] 
g,e-l 
pCt) = pate se e o Bt (AZ 4) 
Ce) 


where ['(g) is the gamma function of g. This expression is the 
probability density function for a gamma distribution of order g with 


a scaling factor 6. 


If x, belongs to a population of uniformly distributed 
a 


integers, the values of y, are limited to 
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and the values of Ly must be derived from these Yas This means that 
the value of any interpulse interval created from a series of uniformly 
distributed integers in the above manner must be one of a specific set 
of values. No other interpulse intervals are permitted. In addition, 
interpulse intervals which are infinite or equal to zero must be 
eliminated if the pulse erate is to be physically realizable. 

The see spectra for a set of such pulse trains with values 
of g ranging from 1 to 8 were estimated from modified periodograms 
according to the methods outlined in Chapter 3 and compared with the 
corresponding theoretical spectra calculated from Equation (4,2), 


Differences between the actual and theoretical spectra were negligible. 
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APPENDIX 3 
FLOW CHARTS FOR COMPUTER PROGRAMS 


The following terms, used in the flow charts, are entered on 


the teletype when required: 
Data Cobtecbion ftogram (Figure 43.1) 

ITIME Sampling interval (msec) 

NPTS Number of points per data block 

NRECS Number of data blocks to be collected 

Sampling Control Commands 

AB Abort 

RS Restart 

GO Begin sampling 

MI’ Stop 

PA Pause--temporarily stop sampling until CO command 

received. Maintain scope display. 

CO Continue 

Data Analysis Program (Figure A3.2) 

RECORDS TO DISPLAY Identification numbers of first and 
last members of a string of consecu- 
tive data blocks. Transfer function 
estimates are calculated by averaging 
spectral estimates based on these 
consecutive data blocks. 

LABEL Two numbers which will appear in the lower left 

corner of the plot to aid in identification. 
Display Commands--Indicate which plot is to be displayed on 


scope. 
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Data Storage Program (Figure A3.3) 


LABEL Any 72 character identifying label to be stored 


Gain 
Phase 
Coherence 
New Plot 
Stop 


Restart 


on the magnetic tape together with the data 


NRECS Number of data blocks to be stored 


Curve Fitting Program (Figure A3.4) 


RECORDS TO BE DISPLAYED 


Display Commands 


RS Restart 


PH Phase plot 


NP New plot 


SL otop 


See data analysis program 
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A3.1 Flow Chart for Data Collection Program 
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APPENDIX 4 
STIMULUS AND SAMPLING CONTROL CIRCULTRY 

The circuit in Figure A4.1 provides bursts of pulses to 
trigger a Hewlett-Packard 5610A analog to digital converter so that it 
collects samples in groups of 128. In addition, gate and trigger sig- 
nals for an Ortec 4710 stimulator are provided so that stimulation 
occurs only during these bursts. 

The sine wave appearing at the base of QQ is converted to a 
pulse train by means of the shaping circuit consisting of QQ and Q,- 
These pulses are coupled through G, and G, to the base of Qe and to the 


it 2 


imput of the counting icircuit mhich consists of Cc and the four J-K 


is 


4 


closes and counting stops. 


flip flops. This situation continues as long as the Q output of FF 
high. When this output changes state, CG, 
When a reset pulse occurs at the base of Qs the counter and flip flops 
are reset, Q becomes high, and pulses may once again pass through G,- 


When 128 pulses have occurred, G, is closed, and counting once again 


2 
ceases. Thus, pulses appearing at the emitter of Qe will occur in 
bursts of 128. The pulse frequency is equal to the frequency of the 
sine wave, and the interval between bursts is equal to the interval 
between reset pulses. G, insures that output pulses do not occur during 


1 


the reset period. 


Trigger pulses for the stimulator are derived from pulses 
applied to the base of Q- These pulses are reshaped by the circuit 
consisting of Q, and Qy> and applied to one input of G, which is open 
whenever the Q output of FF) is high. In this: way, external trigger 
pulses for the stimulator are gated by CG. If the stimulator is to be 


triggered internally it may be gated directly with the Q output of FF): 
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APPENDIX 5 
THE CURVE FITTING TECHNIQUE 
P ; 
If H is a transfer function estimate determined at frequencies 
- 


f,,---f, and H is a function of one or more parameters which is to 


be fitted to this estimate, the parameter values for H which result in 


the best fit are determined when the quantity ¢ is minimized. 


INF 
€= 5 lee j-ucge (2 A5.1 
ok Jf D-HGE,) | (A5.1) 


For single motor units, H is of the form 


kf 
aes = n oe 2TET (A5..2) 


ae 
eee eden 


€ is therefore a function of four parameters Tt, k, © and fo: 

The parameter values for which € is minimized may be 
calculated with subroutine FMFP from the IBM scientific subroutine 
package [136]. This subroutine, which uses the Fletcher-Powell method 
[110] for minimizing a function, requires a user-written subroutine 
which calculates the function value and gradient vector for a given set 
of parameter values. The function value may be determined by 
evaluating Equation (A5.2) at the desired frequencies and substituting 
the results into Equation (A5.1), but evaluation of the gradient vector 
requires a set of expressions for the partial derivatives of € with 


respect to each of the four parameters. 


”~ 
nz” — 


H and H are complex numbers with real and imaginary parts 


Boe a. and Hy Hy respectively. For any parameter \v, 


: é % aH, 7 er 
eas Do _+/7y _Y INS). 
Cea lay H,) Dee (H_-H_)] 
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H can be written as 


H = AC-BD+j (BC+AD) (Aga) 
kE* (£°-£7) 
where A = 7 or (A5.5) 
(£°-F*) "44076767 
n n 
3 
=o E 
re el 0 ee (A5.6) 
(f°-£°) 440-4 £2 
n n 
C = cos 2T£T (A5.7) 
D = -sin 27TfT (A5.8) 
The partial derivatives are 
oH 
Reg eC ok 5 OD Ly OE 
ae A oer sie are B ree D Sai (A5.9) 
dH 
Dee ee 30 ga es 
Oe ee iy a eas Met) 


Expressions for the partial derivatives of A, B, C, and D 
with respect to each of the four parameters T, k, G and a are pre- 
sented in Table A5.1. In order to determine the gradient vector, these 
expressions must be evaluated at the desired frequencies. The results 
may be used in Equations (A5.9) and (A5.10) and a value for each 
element of the gradient vector may then be determined from Equation 
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APPENDIX 6 


DIGITAL REALIZATION OF SECOND ORDER LOW-PASS FILTER 

The frequency response of an analog filter with transfer 
function H(s) is found by evaluating H(s) at points along the imaginary 
axis of the s plane. The equivalent operation for a digital filter 
with transfer function H(z) consists of evaluating this function at 
points along the unit circle in the z plane where z = aoe and T is the 
sampling interval. The imaginary axis on the s plane may be mapped onto 
the unit circle on the z plane by means of the bilinear transformation 


LST} 


z-1 
bs) rey (A6.1) 


If the analog frequency variable is defined as Wy and the 


digital frequency variable is wT, the functions H(w,) and H" (wp T) 


D 


take on the same values when [112] 


D 
= — A6,2 
w, = tan 5 (A6.2) 


Therefore, a second order low-pass digital filter with natural 
frequency wl is equivalent to an analog filter with a natural frequency 


of 
(MAL 
Ws = Thanks (A6 .3) 
@ 2 


Using Equations (A6.1) and (6.2) the transfer function for the digital 


Fagiter, 1s 
iz, 
nals 
H"(z)"= (A6.4) 
z-1.2 z-1 f) 
ee W (——)4+W 
CrP H25 « GFP c 
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This is equivalent to 


it 2 


w (1+2z “+z “) 
H'(z) = Cray ee a oan Oe ee (A6.5) 


2 -1 Zz -2 2 
(1+2zw tw) +z (2u (72) +2 (1-27u tw.) 


and this expression is equal to the ratio of the z transforms of the 
output and input signals y(nT) and x(nT) where n=0,1,2... 


Y 
Bu(ay 2 we (A6 .6) 


If this relationship is substituted into Equation (A6.5) the result is 


Y(z) = 7 [¥(2) (-Bz Cz ")+w2X(2) (14227427) | (A6.7) 
where ee 14200 +0 (A6.8) 
ra 2 
B = 2-1) 


C.= 1-200 tu. 


The inverse z transform of Equation (A6.7) is 


iy (ual) tT) — Cyl m2) Tit + 


y(nT) = = 


(A6.11) 


w 2 


ay {x[nt]+2x[ (n-1)T]+<[ (n-2)T]} 


With this equation, the value of the filtered signal at sampling 
instant n may be calculated when the values of the input signal and the 
filtered signal for the previous two sampling instants are known. Also 
required is the value of the input signal at sampling instant n. When 
“n equals 0 or 1, the previous two sample values are not available, and 
a common practise is to set the unknown values equal to zero. This is 
equivalent to multiplying the continuous signals x(t) and y(t) by unit 


step functions before they are sampled. 
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APPENDIX 7 
THE HUXLEY-JULIAN MODEL 

The following is a summary of the Huxley-Julian model for 
muscle activation and contraction as presented by Stein and Wong [124]. 

Huxley [114] proposed that muscle contraction is produced 
when cross-bridges on the A filaments link with specific sites on the 
I filaments. Because cross-bridges can only become linked on one side 
of their equilibrium position, force is generated in only one 
direction. The time course of interacting cross bridges was described 


by a first-order differential equation. 


ci eS 7 
a (1=-n) f-ng CAT. 1) 


where £ and g are position-dependent variables for making and breaking 
cross-bridges respectively, and n is the proportion of cross-bridges 
that are made. 

The force generated by all the interacting cross-bridges at 


any time, t, is a where 
CO 


P = fk nu du (A7.2) 
& ue ¥ 
Ko is a normalization constant, and u is a normalized variable 
representing the distance from the equilibrium position. The infinite 
limits on the integral imply that integration is extended in both 
directions from the equilibrium position until no longer changes 


appreciably. The equivalent stiffness of the interacting cross-bridges 


is k where 
& 


k = Jk n du (A723) 
g 6) 


Sr Oo. 


The variable f in Equation (A7.1) was assumed to depend on 
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